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I.  INTRODUCTION 

This  report  summarizes  work  carried  out  during  the  past  two 
years  on  two  problems  arising  from  the  interactions  between  acoustic 
fields  and  non-uniform  flow  fields.  Both  problems  are  associated  with 
the  more  general  subject  of  the  stability  of  acoustic  waves  in  solid 
propellant  rocket  motors.  One  purpose  of  the  present  work  was  to 
demonstrate  what  could  be  learned  with  experiments  performed  at  room 
temperature.  It  is  important  to  distinguish  those  problems  which  are 
purely  fluid  mechanical  and  therefore  which  can  be  studied  without  the 
serious  complications  accompanying  combustion  processes. 

In  Section  II  we  describe  the  first  measurements  of  the  loss  of 
acoustic  energy  accompanying  the  flow  of  gas  inward  through  the  lateral 
boundary  of  a  chamber  containing  acoustic  waves.  The  flow  entering 
normal  to  the  surface  must  turn  to  the  axial  direction  parallel 
to  the  boundary;  simultaneously  each  fluid  element  of  the 
mean  flow  must  acquire  acoustic  energy.  It  is  this  second  process, 
involving  action  of  the  unsteady  field  in  the  chamber  upon  the  entering 
flow,  which  appears  as  a  loss  of  energy  for  the  existing  acoustic  field. 

Consider,  for  example,  the  case  of  a  longitudinal  acoustic  field 
in  a  closed  chamber.  For  our  purposes  we  may  assume  that  there  are 
no  viscous  or  other  losses,  so  the  field  will  exist  forever  if  nothing 
else  happens.  But  now  imagine  that  some  gas  is  leaked  in  at  the  lateral 
boundary.  If  there  is  no  way  for  the  gas  to  escape  the  chamber,  the 
average  pressure  and  density  rise  and  the  amplitude  of  the  acoustic 
field  must  decrease  in  order  to  conserve  acoustic  energy.  Let  m., 
be  the  initial  values  of  total  mass  and  density  of  acoustic  energy 


Classically  one 


in  the  chamber,  and  let  m^,  be  the  final  values, 

obviously  has,  for  the  conservative  system, 

m  .3.  =  m,<2, 

11  f  f 

and  because  >  rm,  therefore  3^  <  3^.  But  3 ^  is  proportional  to 

the  square  of  the  amplitude  of  the  pressure  oscillation,  so  the  amplitude 
is  reduced  by  the  addition  of  mass. 

The  same  argument  works  in  reverse:  if  the  bit  of  mass  is 
removed,  the  amplitude  will  be  increased  to  its  initial  value.  Now  in 
steady  flow,  each  element  of  mass  entering  is  matched  by  an  equal 
amount  departing.  Thus  the  preceding  argument  suggests  that  a  steady 
mean  flow  should  not  disturb  the  acoustic  field.  The  conclusion  is  in¬ 
correct  and  the  argument  fails,  because  in  the  real  case  the  influences 
of  viscous  forces  cannot  be  ignored.  For  a  non-uniform  flow  field,  the 
forces  required  to  accelerate  the  incoming  fluid  necessarily  involve  in¬ 
elastic  processes  and  energy  is  dissipated.  In  fact,  half  of  the  work 
done  by  the  gas  in  the  chamber  on  the  incoming  fluid  is  dissipated  as 
heat  and  half  provides  the  acoustic  energy  acquired  by  the  entering  mass 
(Ref.  1).  As  a  result  there  is  necessarily  a  loss  of  energy  associated 
with  the  inward  flow  of  gas  which  initially  has  no  acoustic  energy. 

That  energy  loss  -  a  loss  for  the  acoustic  field  existing  in  the 
chamber  -  is  often  referred  to  as  the  'flow  turning  loss'  because  turning 
of  the  average  flow  is  necessarily  involved.  It  is  really  a  rate  of  loss, 
proportional  to  the  velocity  of  the  incoming  flow  perpendicular  to  the 

1.  Culick,  F.  E.  C.  "Remarks  on  Entropy  Production  in  the  One- 
Dimensional  Approximation  to  Unsteady  Flow  in  Combustion 
Chambers",  Comb.  5ci.  and  Tech.  .  V.  15  (1977),  pp.  93-97. 


boundary.  Its  existence  was  first  predicted  theoretically  (Ref.  2)  in 
computations  based  on  the  one-dimensional  approximation  but  has  not 
previously  been  measured.  The  measurement  technique  used  in  the  present 
work  is  an  extension  of  that  used  in  earlier  work  (Refs.  3,  4  and  5). 

A  resonance  tube  with  provision  for  introducing  and  exhausting  an 
average  flow  is  the  central  apparatus  (see  Figure  1,  p.  19).  Flow  enters 
through  the  porous  faces  of  pistons  fitted  to  the  ends  of  the  tube.  The  pistons 
are  driven  to  excite  standing  acoustic  waves.  The  flow  exhausts  through  a  vent 
in  the  center  of  the  tube,  enforcing  symmetry  of  the  average  flow  field.  For 
measurements  of  the  flow  turning  losses,  cylindrical  porous  tubes  are  mounted 
at  intermediate  positions  between  the  ends  and  the  center.  In  principle, 
the  procedure  is  straightforward.  The  forcing  frequency  is  varied  across 
the  fundamental  resonant  frequency  of  the  system,  and  the  slope  of  the 
resonance  curve  is  determined.  The  width  of  the  curve  is  proportional 
to  the  net  losses  in  the  system. 

In  the  previous  work,  the  losses  associated  with  the  porous  pistons 
and  the  gains  associated  with  the  exhaust  vent  were  determined  over  a 

2.  Culick,  F.  E.  C.  "The  Stability  of  One-Dimensional  Motions  in  a 
Rocket  Motor",  Comb,  and  Sci,  and  Tech.  ,  V.  7  (1973),  pp. 

165-175. 

3.  Magiawala,  K.  "Measurements  of  Energy  Exchange  Between  Acoustic 
Fields  and  Non-Uniform  Flow  Fields",  Ph.  D.  Thesis,  California 
Institute  of  Technology  (May  1978). 

4.  Culick,  F.  E.  C.  and  Magiawala,  K.  ,  "Measurements  of  Energy 
Losses  Associated  with  Interactions  Between  Acoustic  Waves  and  a 
Steady  Flow  Field",  Air  Force  Rocket  Propulsion  Laboratory, 

Report  AFRPL-TR-78-6  (March  1978). 

5.  Magiawala,  K.  and  Culick,  F.  E.  C.  "Measurements  of  Energy 
Exchange  Between  Acoustic  Fields  and  Non-Uniform  Steady  Flow 
Fields",  J.  Sound  and  Vib.  (to  be  published). 
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range  of  average  Mach  number.  Those  results  are  used  here,  so  the 
flow  turning  loss  is  then  found  by  performing  additional  tests  with  the 
cylindrical  elements  in  place. 

Serious  difficulties  arise  because  the  quantities  to  be  measured 
are  relatively  small  and  are  determined  indirectly  by  forming  the  dif¬ 
ferences  of  relatively  large  numbers.  Thus  great  care  must  be  taken 
to  achieve  as  high  precision  as  possible.  The  procedures  followed  have 
been  described  fully  in  References  3  and  4.  We  encountered  further 
obstacles  here  in  determining  the  admittance  functions  for  the  cylindrical 
samples.  There  are  two  reasons.  First,  we  have  been  unable  to  devise 
a  means  of  measuring  the  admittance  of  a  curved  sample,  so  we  have 
been  forced  to  assume  that  the  values  are  the  same  as  those  for  a  flat 
sample.  Second,  it  is  a  consequence  of  the  way  in  which  the  materials 
are  manufactured  that  we  have  been  unable  to  obtain  flat  samples  having 
the  same  porosity  as  the  cylindrical  samples.  We  have  therefore  used 
values  extrapolated  from  those  actually  measured.  The  technique  is 
described  in  Section  2.  5. 

Those  problems  have  contributed  substantially  to  experimental  un¬ 
certainties  in  the  final  results.  Our  measurements  have  shown  that  the 
flow  turning  loss  does  indeed  increase  with  the  average  flow  speed,  but 
is  approximately  three  times  the  value  predicted  by  the  approximate  one¬ 
dimensional  analysis.  The  experimental  errors  are  quite  large.  We 
have  reduced  the  measurement  errors  as  far  as  possible,  and  it  appears 
that  the  only  way  to  reduce  the  uncertainties  in  the  final  results  is  to 
use  larger  apparatus  with  higher  flow  rates.  The  absolute  numerical 
values  to  be  measured  would  therefore  be  larger,  thereby  reducing  the 


relative  sizes  of  the  experimental  errors.  This  conclusion  applies  to 
all  experimental  work  with  measurements  for  problems  involving  inter¬ 
actions  between  acoustic  waves  and  non-uniform  flow  fields. 

In  Section  III  we  describe  measurements  of  the  loss  of  acoustic 
energy  due  to  the  action  of  a  choked  exhaust  nozzle.  This  influence, 
represented  by  an  admittance  function,  is  a  substantial  contribution  to 
attenuation  of  longitudinal  waves  in  a  rocket  motor.  The  influence  of  a 
de  Laval  nozzle  was  first  examined  theoretically  long  ago  (Refs.  6 
and  7).  The  first  measurements,  using  apparatus  with  steady  flow  (Ref.  8) 
verified  theoretical  predictions  but  were  limited  to  simple  nozzles. 

The  admittance  function  was  determined  by  direct  measurement  of  the 
velocity  fluctuations  with  a  hot  wire  anemometer,  and  the  pressure 
oscillations  with  a  static  pressure  probe.  By  far  the  most  extensive 
measurements  were  those  reported  in  References  9-12  using  a  blow-down 
facility.  Some  results  were  obtained  (Ref.  11)  for  the  admittance 
of  a  nozzle  with  secondary  flow. 

6.  Tsien,  H.  S.  "The  Transfer  Function  of  Rocket  Nozzles",  ARSJ, 

V.  22,  No.  31  (May- June  1952),  p.  139. 

7.  Crocco,  L.  and  Cheng,  S.  -T.  Theory  of  Combustion  Instability  in 
Liquid  Propellant  Rocket  Motors,  AGARDograph  No.  8,  Butterworths 
Publications,  Ltd.,  London  (1958)  Appendix  C. 

8.  Crocco,  L.  ,  Monti,  R.  and  Grey,  J.  "Verification  of  Nozzle  Ad¬ 
mittance  Theory  by  Direct  Measurement  of  the  Admittance  Parameter", 
ARSJ.  V.  31,  No.  6  (June  1961)  pp.  771-775. 

9.  Zinn,  B.T.  et  al  "Experimental  Determination  of  Three-Dimensional 
Liquid  Rocket  Nozzle  Admittances",  AIAA  J. ,  V.  11,  No.  31 

(Mar.  1973),  pp.  267-272. 

10.  Bell,  W.  A.,  Daniel,  B.  R.  and  Zinn,  B.T.  "Experimental  and 
Theoretical  Determination  of  the  Admittances  of  a  Family  of 
Nozzles  Subjected  to  Axial  Instabilities".  J.  Sound  and  Vib. 

V.  30,  No.  2  (1973),  pp.  179-190.  ’ 
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Although  the  method  used  in  Ref.  8  seems  to  have  given  reasonable 
results,  there  are  considerable  difficulties  associated  in  measuring  the 
small  acoustic  velocities,  including  distortion  and  low  signal/noise  ratio. 
Moreover  many  measurements  should  be  made  over  a  plane  near  the 
entrance  to  the  nozzle,  and  averaged,  in  order  to  give  a  proper  represen¬ 
tation  of  the  influence  of  the  entire  nozzle.  That  procedure  was  not 
followed  in  Ref.  8.  This  technique  would  likely  become  particularly 
tedious  and  time  consuming  for  a  nozzle  with  secondary  flow.  In  that 
case  considerable  effort  would  likely  be  required  to  determine  the  best 
reference  plane  for  defining  the  admittance  function. 

The  data  reported  in  Refs.  8-12  seem  to  be  quite  ieasonable,  although 
there  is  no  discussion  of  experimental  errors  or  the  cost  of  tests.  It  appears 
that  the  cost  of  using  a  blow-down  facility  for  this  purpose  might  be 
quite  high.  The  section  upstream  of  the  nozzle  is  used  as  an  impedance 
tube,  so  the  acoustic  pressure  must  be  measured  at  many  locations. 

For  transient  tests,  many  pressure  transducers  are  required,  and  data 
processing  becomes  a  major  expense  because  of  the  equipment  required. 

A  primary  purpose  of  the  present  program  has  been  to  determine 
how  well  one  might  use  a  small  continuously  operating  supersonic  wind 
tunnel  to  measure  the  nozzle  admittance.  The  main  advantages  are  the 
cost  ($25  per  hour  for  the  tunnel)  and  the  flexibility  of  continuous  operation 

11.  Zinn,  B.T.  et  al  "Damping  of  Axial  Instabilities  by  the 
Minuteman  II,  Stage  III  and  Minuteman  III,  Stage  III,  Exhaust 
Nozzles",  Georgia  Institute  of  Technology,  Report  AFRPL-TR- 
72-71  (August  1972). 

12.  _ _  "Nozzle  Design  Considerations  for  Attenuation 

oi  Axial  Instabilities  in  the  Minuteman  II  and  III,  Stage  III 
Rocket  Motors",  Georgia  Institute  of  Technology,  Report 
AFRPL-TR- 73-69  (September  1973). 
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which  permits  studying  various  aspects  of  the  problem  with  relative  ease. 

For  example,  the  spatial  characteristics  of  the  average  and  unsteady 
flow  field  are  particularly  important  because  to  reduce  the  data,  as  in 
Refs.  9-12  and  here,  the  assumption  is  used  that  the  flow  fields  are 
uniform  over  a  cross-section  of  the  tube.  The  mean  flow  obviously 
cannot  be,  because  of  the  growth  of  the  boundary  layers.  How  significantly 
this  approximation  affects  the  final  results  has  yet  to  be  established;  time 
did  not  permit  examining  the  problem  in  this  work.  However,  we  have 
obtained  some  results  showing  that  the  acoustic  field  is  not  one-dimensional 
under  all  conditions. 

It  is  also  possible  with  the  techniques  used  here  to  obtain  the  necessary 
measurements  with  only  two  pressure  transducers  and  therefore  only  two 
channels  of  data.  Overall,  then,  the  cost  of  the  equipment  is  much  below  that  of 
a  blow-down  facility.  The  quality  of  the  data  should  be  better  -  i.  e.  one  should 
be  able  to  reduce  the  experimental  errors  because  the  techniques  are 
simpler  with  steady  operation  -  but  this  has  not  been  proved  with  the 
present  work. 

In  principle,  the  use  of  the  subsonic  section  of  a  wind  tunnel  as 
an  impedance  tube  is  really  quite  straightforward.  We  had  assumed,  on 
the  basis  of  the  work  reported  in  Refs.  9-12,  that  the  well-known  technique 
based  on  measurements  of  the  pressure  amplitude  would  be  directly 
applicable.  After  considerable  effort  and  time,  we  concluded  that  the 
assumption  is  wrong.  It  happens  that  if  the  losses  at  the  end  of  the 
impedance  tube  are  relatively  small,  as  the  case  is  for  choked  nozzles, 
it  is  essential  that  both  the  amplitude  and  phase  of  the  acoustic  pressure 
be  measured.  Both  pieces  of  information  must  be  used  to  obtain  accurate 
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results  for  the  real  and  imaginary  parts  of  the  admittance  function. 

This  problem  is  discussed  in  Sections  3.2  and  3.6. 

Because  we  have  examined  the  influences  of  secondary  flow  and 
submergence,  the  configurations  tested  have  shapes  for  which  the  admit¬ 
tances  cannot  be  calculated  with  any  available  analysis.  The  values 
inferred  from  our  measurements  seem  to  be  reasonable,  with  the 
possible  exception  of  an  apparent  energy  gain  at  the  highest  frequency 
used  (2000  Hz).  That  may  be  due  to  resonance  within  the  plenum 
chamber  surrounding  the  nozzle,  a  provision  required  for  the  introduction 
of  secondary  flow.  In  any  case  the  experimental  uncertainty  is  large. 

The  signal/noise  ratio  causes  problems  which  increase  with  the 
Mach  number  of  the  average  flow.  Up  to  M  =  0.  04  a  tracking  filter 
(bandwidth  =  0.  5  Hertz)  is  entirely  adequate.  The  values  for  the  ad¬ 
mittance  function  can  be  determined  with  less  than  5#  error.  For 
M  >  0.  1  the  problem  is  more  serious  in  the  vicinity  of  pressure  minima. 
The  signal/noise  ratio  is  of  course  much  reduced,  but  results  could  be 
extracted  by  numerical  averaging  (which  we  have  not  done  in  this  program). 
A  more  serious  difficulty  we  encountered  was  apparent  fluctuation  of  the 
axial  location  of  the  pressure  minimum.  Time  did  not  permit  thorough 
examination  of  the  problem,  and  we  cannot  state  whether  it  is  soluble 
simply  by  averaging  the  data.  If  it  cannot  be  handled  in  that  manner, 
then  the  use  of  the  tunnel  as  an  impedance  tube  may  be  blocked  for 
higher  Mach  numbers. 

Should  that  be  the  case,  then  one  would  be  forced  to  return  to 
direct  measurements,  as  in  Ref.  8.  However,  the  best  technique  for 
measuring  the  velocity  is  now  laser-Doppler  anemometry.  There  are 


certain  practical  problems  when  circular  transparent  test  sections  must 
be  used,  but  it  appears  that  this  might  be  a  fruitful  course  to  follow 
in  any  case.  The  importance  of  nozzle  losses  in  the  problem  of  com¬ 
bustion  instability  demands  that  a  method  be  developed  which  is  relatively 
cheap;  reliable;  and  can  be  used  routinely  to  measure  the  acoustical 
properties  of  nozzles  used  in  operational  rocket  motors. 
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H.  MEASUREMENTS  OF  FLOW  TURNING  LOSSES 

Since  the  existence  of  the  flow  turning  loss  was  first  predicted  in 
Ref.  1,  there  has  been  a  certain  amount  of  controversy  concerning  both 
the  qualitative  and  quantitative  truth  of  the  result;  see  Ref.  13,  for 
example,  for  a  brief  discussion.  The  data  reported  here  constitute 
the  first  evidence,  obtained  from  purely  acoustics  experiments,  showing 
that  the  loss  does  arise  in  real  flows.  It  is  therefore  appropriate  to 
begin  here  with  a  short  resume  of  the  theoretical  origins,  after  which 
we  describe  the  experimental  work  carried  out  in  the  present  program. 

2.  1  Analysis  of  the  F low  Turning  Loss 

The  flow  turning  loss  first  appeared  in  the  analysis  (Ref.  1)  of 
longitudinal  acoustic  waves  in  a  combustion,  using  the  one -dimensional 
approximation.  This  is  essentially  an  extension  to  unsteady  flow  with 
combustion  of  the  problem  of  flow  in  a  duct  with  mass  addition  (e.  g. 
Ref.  14,  Chapter  8).  The  pressure  loss  in  flow  with  mass  addition 
(see  Table  8.  1  of  Ref.  14)  corresponds  to  the  loss  of  acoustic  energy 
associated  with  mass  addition  in  the  present  case. 

We  can  demonstrate  the  result  by  considering  the  simplest 
problem  of  flow  in  a  uniform  channel  with  no  combustion.  Within 
the  one-dimensional  approximation,  the  general  conservation  equations 
and  the  equation  of  state  can  be  written 


13.  Culick,  F.  E.  C.  (Ed.),  "T-Burner  Testing  of  Metallized  Propellants", 
AFRPL  Report  TR-74-28  (Oct.  1974). 

14.  Shapiro,  A.  H.  The  Dynamics  and  Thermodynamics  of  Compressible 
Fluid  Flow.  The  Ronald  Press  Company,  New  York  (1953). 
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f?  +  h  (pu)  =  w  <2- 

jt  <pu>  +  -h  (pn2)  +  If =  usw  <2-2) 

(peo)  +  ~k  (pueo)  +  ~k  (pu)  =  hosw  (2-3) 

p  =  pRT  (2.4) 


where 

w  =  ~  /  mbdq  (2.  5) 

°c 

The  source  term  on  the  right  hand  side  of  (2.2)  represents  the  initial 

axial  momentum  of  the  gas  flowing  in  at  the  boundary,  and  the  right 

hand  side  of  (2.3)  represents  the  enthalpy  (internal  energy  plus  p-v 

work)  carried  in.  The  mass  addition  at  the  boundary  is  m  (mass/area- 

second);  the  perimeter  is  q;  and  Sc  is  the  cross-sectional  area. 

2 

Subscript  zero,  (  )^,  denotes  the  stagnation  value:  e^  =  e  +  u  /2  = 

2  - 
CyT  +  u  /2  and  hQ  =  e  +  p/p  +  u2/2. 

Equation  (2.2)  can  be  re-written,  with  use  of  (2.1): 


9u  9u  9p  .  . 

p  di  +  pu  97  +  =  (us  "  u)' 


(2.6) 


The  right  hand  side  of  this  equation  represents  the  force  required 
to  accelerate  the  added  mass  to  the  local  axial  speed  u.  That  this  is 
an  inelastic  and  therefore  dissipative  process  may  be  understood  in  an 
elementary  way  from  the  equations  of  mechanics.  A  clear  and  extended 
discussion  of  the  point  may  be  found  in  Ref.  15. 


15.  Rauscher,  M.  Introduction  to  Aeronautical  Dynamics,  John  Wiley 
and  Sons,  Inc.,  New  York  (1953)  p.  50. 


Similarly,  (2.3)  can  be  re-formed  and  the  kinetic  energy  removed 
by  subtracting  u  times  the  last  equation  to  give 

9T  /-  9T  9u  tv.  \  i  \  r>  n\ 

P°v  ~dt  +  PuCv  aT  +  p  to  =  (h0s  -  eo)w  -  U(US  -  u)w  (2-  7) 

We  have  chosen,  as  in  earlier  works,  to  treat  the  mass  addition 
and  its  influences  as  sources  averaged  over  the  cross-section  of  the 
channel.  One  could  also  begin  with  the  equations  for  two-  or  three- 
dimensional  flow  and  average  over  the  cross-section.  The  influences  of 
the  sources  at  the  boundary  then  arise  through  the  boundary  conditions. 
With  some  assumptions  concerning  definitions  of  averaged  quantities 
("closure  conditions")  the  same  equations  (2.  l)-(2.3)  can  be  recovered; 
some  details  of  this  procedure  may  be  found  in  Ref.  16. 

It  is  convenient  to  use  the  equation  for  the  pressure,  constructed 
by  adding  Cv  times  the  continuity  equation  (2.1),  to  (2.7)  and  using 
the  equation  of  state  (2.4);  the  final  result  is: 

-|^  +  yp-|^  +  u-|^  =  [a2  +yRAT  +-|-(u2-u2)]  w  +  u(u  -  ug)w  (2.8) 

where  AT  =  Tg  -  T  is  the  local  temperature  difference  between  the 
bulk  flow  in  the  channel  and  the  entering  flow. 

We  now  assume  that  the  flow  enters  with  no  axial  motion,  ug  =  0, 
and  with  temperature  equal  to  the  local  temperature  in  the  duct,  so 
Tg  =  T;  these  assumptions  apply  to  the  experiments  described  below. 
Equations  (2.  6)  and  (2.  8)  become 

16.  Van  Moorhem,  W.  K.  "An  Investigation  of  the  Origin  of  the  Flow- 

Turning  Effect  in  Combustion  Instability",  17th  JANNAF  Combustion 
Meeting  (September  1980). 
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P  dt  +  pu 


9u 


dz 


dz 


(2.9) 


,  9u  3p  2  R  1.  2 

+  *>  al  +  u  at  =  a  w  +  (C~  •  2)u  w  (2-  10) 

V 

Analysis  of  the  stability  of  small  disturbances  is  based  on  the  linearized 

forms  of  these  two  equations;  set  u  =  u  +  u'  etc.  ,  and  ignore  terms 
2  —2 

of  order  u'  and  u  to  find 


9 £ 
at 


p 


9u' 

at 


(uw) ' 


(uu1) 


(2. 11) 


-  9u' 


yp 


dz 


=  (aw)' 


yp' 


du 

d2 


(2.  12) 


We  assume  that  while  the  average  speed  u  varies  with  position,  the 
average  density,  temperature  and  therefore  pressure  are  constant. 

Now  form  the  wave  equation  for  the  pressure  fluctuation  by  differ¬ 
entiating  (2.12)  with  respect  to  time  and  substituting  (2.11);  witn  some 
rearrangement  we  find 


8  » i 

al(=2  w) 

cl 


_  g2  _  y_  9p'  du 

p  — j  (uu')  +  ,  at  dz 

az  a^ 


+ 


ataz 


(2.  13) 


For  our  purposes,  to  demonstrate  the  appearance  of  the  flow  turning 
loss,  only  the  first  terms  on  the  right  hand  sides  of  (2.  11)  and  (2.  13) 
are  required.  The  wave  equation  (2.  13)  with  the  boundary  condition  to 
be  applied  at  the  ends  of  the  channel,  given  by  (2.  11)  are 


_a_ 

az 


(uw)' 


(2.  14) 


-^7  =  -  (uw)'  (z=0,  L) 


(2.  15) 


An  arbitrary  disturbance  treated  within  linear  analysis  can  be 

represented  as  a  Fourier  synthesis  of  sinusoidal  waves.  Stability  is 

assured  only  if  all  contributions  over  the  entire  frequency  range  are 

stable.  Hence  we  can  treat  stability  by  examining  the  typical  wave 

having  frequency  uu ,  and  write 

iakt 

p'  =  p  e  (2.  16) 

where  k  =  (uu  -  ia)/a  is  the  complex  wavenumber.  Substitute  (2.  16) 
into  (2.  14)  and  (2.  15)  to  find 

+  k2  =  -  —  (uw)  (2.  17) 

dz2  dz 

ff  =  -  (uW)  (z=0,  L)  (2.18) 

We  shall  give  uw  explicitly  below. 

For  stability,  the  imaginary  part,  a,  must  be  negative.  This 
represents  a  loss  of  acoustic  energy.  It  is  a  property  of  this  linear 
system  that  if  the  time  averaged  total  energy  is  (<$) ,  proportional 
to  the  square  of  the  pressure  fluctuation,  then  the  constant  a  is 
related  to  the  rate  of  change  of  energy  by 

2a  =  w  ^dp  (2- 19) 

Hence  the  problem  comes  down  to  calculating  the  decay  or  growth 
constant,  a,  due  to  the  perturbations  represented  by  the  right  hand  sides 
of  (2.  14)  and  (2.  15).  This  is  done  most  directly  by  comparing  the  problem 
containing  perturbations  with  the  unperturbed  problem.  For  the  latter 
case,  the  mode  shape  p^  and  wavenumber  k^  are  governed  by  the 
equations 
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,2- 

P*  *  - 
- T  +  k«  Pi 


dz 


± l 
dz 


=  0 


=  0  (z=0,  L) 


(2.  20) 


(2.  21) 


Multiply  (2.  17)  by  p^,  (2.20)  by  p,  subtract  the  two  equations,  and 
integrate  over  the  length  of  the  chamber,  0  <  z  <  L.  to  find 

L  L 


^  2  ^  ^  * 

/  Cp  ^  ^-|  -  P  - =r]dz  +  (k2-k^)  /  pp-dz  =  -  /  JT  (uw)  P  .dz 

0  *  dz  dz'1  Z  0  1  0  d  £ 


Now  integrate  the  first  term  by  parts  and  use  the  boundary  conditions 


(2.18)  and  (2.21)  to  produce  the  formula  for  k  : 

dz 


2  2  1  \  *. 

k  =  k  +  — j  f  (uw)  -r — 

*  e]  0  dz 


where 


2  k  L>  2 

E.  =  /  P  P,dz  /  p  .dz 

1  0  1  0  1 


(2. 22) 


(2.23) 


Finally,  we  expand  the  product  uw  to  obtain  the  explicit  form 
of  the  flow  turning  loss.  Assume  that  the  mass  flux  inward  is  constant 
over  the  perimeter,  so  w  =  m^q/Sc,  and  linearization  proceeds  as 
follows: 


(uw)1  =  (u  +  u')(w  +  w')  -  uw 

—  uw'  +  u'w  =  (um^  +  u'm^)(q/Sc) 


iakt 

(uw)  e 


Note  that  m^  is  the  fluctuation  of  the  mass  flux  inward  at  the  boundary, 
while  u'  is  the  acoustic  velocity  parallel  to  the  surface.  When  this 
formula  is  substituted  into  the  last  term  of  (2.22),  u'  and  therefore  u 
is,  to  the  order  treated  here,  replaced  by  the  unperturbed  acoustic  velocity. 


i 
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i 


iakt  iakt 

(uw) '  (um  +  u^m^)  -jj-  e  =  (uw)e 

c 

The  formula  (2. 22)  is  therefore 


(2.  24) 


k  =  k 


->  q/S  L  _  L  dp. 

+  -3T  [  L u  Vidz  +  L  Vj  dr dz] 


(2. 25) 


'l 


0 


0 


The  first  term  in  brackets  represents  the  response  of  the  incoming  flow 
to  disturbances  and  will  later  be  written  in  terms  of  the  admittance 
function  for  the  surface.  It  is  the  second  term  which  represents  the 
flow  turning  loss.  For  the  unperturbed  acoustic  field,  classical  acoustics 
gives  the  result  u^  =  (i/p  a  k^)dp^/dz  so  the  integral  can  be  written 

2 

(2.  26) 


dp 

t  —  *  _  i  j 

1  "Vx  d T  dz  = 


L  dp. 

==—  {.  mb(dT>  dz 

P  ak.  0 


Now  k2  =  (uu  -  ia)2  /a2 


(uj/a)2 


i(2aJi)/a),  since  for  small  perturbations, 
the  terms  in  brackets  in  (2.  25)  are  small,  and  therefore  a  «  uu.  With  only  the 
flow  turning  loss  included,  (2.  25)  may  now  be  written 


,w.2  .  ,2auu,  ,2  ,  .  q/Sc 

(_)  .  l(_)  =  k£  +  i 

a 


1 


2 _ 

E .  p  ak 


L  _  dp* 

/  mb(-^-)  dz 


(2.27) 


'l  K  ‘l 

This  formula  can  be  split  into  real  and  imaginary  parts,  the  real  part 
showing  that  the  frequency  is  unchanged  by  the  perturbations,  uu  =  ak., 
and  the  imaginary  part  giving  the  formula  for  the  decay  constant  due 
to  the  flow  turning  loss. 


"ft 


2S  E.  0  p 
cl  y 


L  m,  .  dp .  2 
r  b  .  1  *1. 

1  —  %  — 1 1  d2 


(2.  28) 


It  is  easier  to  remember  this  result  if  it  is  written  in  a  different 

2 

form.  First,  the  constant  is  related  to  the  time-averaged  total 
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acoustic  energy  by  the  following: 


<$>  =  (///  (7  p|u'!2  +  7^^-  )dv> 

P  a 


L  Le:'2 


-  <Sc  {  ~~2 

0  pa 

S  L 


dz  ) 


1  c  7-2 

"  2  -  -2  {  P£dz 

pa  0 

1  S 

1  c  l 

=  i  =-=r 

p  a 


(2.29) 


We  have  used  the  fact  that  the  time-averaged  kinetic  and  potential 
energies  are  equal.  Now  for  the  unperturbed  acoustic  field,  with 
phases  measured  with  respect  to  the  pressure,  so  p^  is  real,  we 


have 


dp« 2  >« 

'dT»  *  'dr 


*  |i  p»ut|2  =  (p  ak^)2 


Thus  the  integral  in  (2.28)  can  be  written,  with  =  p  u^, 

L  _  1  dp  2  L  _  ___  ?  2 

{  Vi^  dT*  dz  =  /  ^  <Pa>  1**1  dz 


After  this  relation  and  (2.29)  have  been  substituted  into  (2.28)  we  have 


2aft  =  -  vr> 


(2.  30) 


where  dS,  =  qdz  is  the  element  of  area  on  the  lateral  boundary, 

b 

The  time-averaged  value  of  the  acoustic  kinetic  energy  density 
is  p|u^|  /4,  the  factor  1/4  composed  of  1/2  from  the  definition 
of  the  kinetic  energy  and  1/2  from  the  time  averaging.  Thus  the 
integral  represents  twice  the  time-averaged  rate  at  which  the  incoming 
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flow  acquires  acoustic  kinetic  energy.  It  is  twice  because  energy  is 
dissipated  (i.  e.  is  associated  with  entropy  production)  at  exactly  the 
same  rate  at  which  the  incoming  flow  gains  its  acoustic  kinetic  energy. 
Hence,  the  acoustic  field  within  the  chamber  suffers  a  fractional  rate 
of  energy  loss,  represented  by  (2.30)  which  is  twice  the  rate  at  which 
mechanical  energy  is  gained  by  the  gas  flowing  in  at  the  boundary.  This 
marvelous  factor  of  two  is  explained  in  Ref.  15  for  a  problem  of  elementary 
mechanics  and  is  verified  by  computation  of  the  rate  at  which  entropy 
is  produced  (Ref.  1  ). 

The  experimental  work  now  to  be  described  is  intended  to  provide 
values  of  which  will  then  be  compared  with  the  values  calculated 

with  the  formula  (2.  28). 

2.  2  Analysis  of  the  Apparatus 

The  resonance  tube  is  sketched  in  Figure  1.  Much  of  the  hard- 
were  is  that  used  in  earlier  work  (Refs.  4  and  5)  to  determine  the 
influence  of  the  exhaust  vent  on  acoustic  waves.  The  new  portions 
are  the  cylindrical  porous  tubes  -  the  "lateral  flow  elements"  -  installed 
as  a  means  of  introducing  flow  in  at  the  lateral  boundary.  Thus,  flow 
in  the  tube  enters  through  both  the  pistons  at  the  ends,  which  serve  to 
excite  acoustic  waves,  and  through  the  cross  flow  elements.  In  Figure  2, 
the  various  flow  and  geometrical  parameters  are  defined. 

In  all  the  tests  reported  here,  the  resonance  tube  was  driven, 
in  its  fundamental  mode,  by  the  oscillating  pistons.  The  acoustic 
pressure  then  has  a  node  at  the  center  and  is  maximum  at  the  ends 
of  the  tube.  When  the  frequency  of  the  driving  mechanism  is  swept 
through  the  resonant  frequency,  the  amplitude  of  the  pressure  oscillation 


i/4M  B8.K  MICROPHONE 


POROUS  TUBE 


Figure  2.  Sketch  Showing  Definitions  for  the  Resonance  Tube 
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describes  a  bell- shaped  curve  with  maximum  very  closely  at  the 
resonant  frequency.  Thus  the  system  appears  to  behave  as  a  simple 
mass /spring/dashpot  system. 

It  is  therefore  possible  to  relate  the  width  of  the  resonance  curve 
to  the  losses  in  the  system.  The  analysis  supporting  that  assertion 
has  been  given  in  Ref.  4  and  need  not  be  repeated  here.  The  main 
result  is  that  if  Auu  is  the  width  of  the  curve  at  the  "half-power 
points",  the  frequencies  at  which  the  amplitude  is  reduced  by  l/Jz" 
from  its  maximum  value,  then  the  net  value  of  the  attenuation  coefficient 
is 

»  =  1  **1/2  <2-31> 


This  formula  is  the  basis  for  the  testing  procedure.  Because  a  is  the 
sum  of  many  contributions,  a  sequence  of  tests  is  necessary  to  determine 
the  separate  parts  until  finally  the  desired  quantity,  in  this  case  a^, 
can  be  determined. 

The  sequence  of  tests  required  is  the  same  as  that  used  in  Refs. 

4  and  5  to  determine  the  influence  of  the  vent,  with  one  additional  test 
to  obtain  the  contribution  from  the  lateral  flow.  From  linear  analysis 

•  jjc 

of  the  configuration  shown  in  Figure  1,  the  decay  constant  a  consists 
of  five  parts. 


a  =  a  +  a.  +  ar,  +  a  +  a , 
p  t  ft  v  d 


(2. 32) 


where  the  terms  on  the  right  hand  side  represent  the  following  con¬ 
tributions: 


a  :  attenuation  constant  associated  with  the  admittance  and  flow  for 
P 

the  faces  of  the  piston  drivers; 

Qp  =  4f(Ajjr)  +  Mp)  (2.33) 


In  the  remainder  of  this  report  we  use  the  sign  convention  that  energy  losses 
are  represented  by  positive  values  of  the  decay  constant  a. 
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Q 


F 


Q 


ft5 


a  : 


v 


a 


d: 


attenuation  constant  associated  with  the  admittance  and  flow 
for  the  cylindrical  lateral  flow  elements; 

\  *  *<Air>  +  st> 


(2.34) 


c 

attenuation  constant  associated  with  the  flow  turning  losses,  to 
be  determined; 

attenuation  constant  associated  with  the  influence  of  the  exhaust 
vent; 

attenuation  constant  associated  with  other  losses,  mainly  radiation 
and  viscous  losses  at  the  lateral  boundary. 


The  formula  (2.33)  is  a  direct  consequence  of  the  one-dimensional 
analysis.  Ref.  1.  In  subsequent  work.  Refs.  13  and  17,  the  result  (2.34)  was 
obtained  in  an  analysis  of  a  variable  area  T-burner;  <3^  is  a  geometrical 
factor  accounting  for  the  finite  extent  of  the  region  through  which  flow 
enters  at  the  lateral  boundary.  For  the  single  configuration  tested  in 
the  present  work,  (3^  =  0.51. 

All  of  the  terms  on  the  right  hand  side  of  (2.32)  are  determined 

experimentally.  The  admittance  functions,  A^1" ^  and  A^1^  are  measured 

by  using  an  impedance  tube,  and  the  Mach  numbers  M  and  M.  are 

P  z 

found  by  measuring  the  flow  rates  through  the  porous  samples.  We  have 
chosen  experimental  conditions  such  that  the  attenuation  constant  due  to 
the  vent  can  be  taken  from  our  earlier  work.  Refs.  4  and  5.  Finally, 
dj  is  determined  as  part  of  the  present  tests  with  the  resonance  tube 
used  to  obtain  a^. 

17.  Culick,  F.  E.  C.  "Linear  Analysis  of  One-Dimensional  Oscillations 
in  a  Variable  Area  T-Burner",  9th  JANNAF  Combustion  Meeting 
(1972). 
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2. 3  Description  of  the  Apparatus  and  Instrumentation 

2.3.  1  Cold  Flow  Resonance  Tube 

Figure  1  shows  the  configuration  of  the  resonance  tube.  The  tube 
itself  is  made  of  brass,  having  inside  diameter  equal  to  3.81  cm,  and 
0.  32  cm  wall  thickness;  for  the  tests  at  445  Hz,  the  length  is  38.  65  cm. 
The  circular  exhaust  vent  is  of  brass,  having  0.32  cm  wall  thickness, 
inside  diameter  equal  to  1.27  cm,  and  is  3.30  cm  long.  It  has  a  plug, 
machined  so  that  when  inserted  the  vent  is  closed  by  a  surface  smoothly 

fitting  the  inside  surface  of  the  main  resonance  tube.  Further  details 

of  the  tube  and  the  dimensions  are  given  in  Figure  2.  The  dimensions 

were  chosen  to  give  p  =  0.22,  PQ  =  0.38,  Sbg/Sc  =  4.53  and  constant 

axial  flow  Mach  number,  for  flow  through  the  porous  pistons, 

M  =  4. 43  X  10'1. 

P 

The  lateral  flow  was  introduced  through  two  porous  tubes  having 
inside  diameter  equal  to  3.81  cm  and  fitted  in  manifolds.  Each  manifold 
has  three  series  of  orifices  which  help  distribute  the  flow  uniformly 
around  the  periphery  of  the  tube. 

The  instrumentation  used  has  been  described  in  Refs.  3-5,  and 
illustrated  in  Figure  3.  The  major  addition  here  is  the  cross-flow 
supply  and  the  related  monitoring  instrumentation  which  is  similar  to 
that  for  the  axial  flow  supply  and  its  monitoring  instrumentation.  To 
avoid  interference  from  airborn  noise  and  local  reflections  within  the 
room,  the  resonance  tube  was  placed  in  an  anechoic  chamber.  The 
air  supply  and  the  instrumentation  were  located  outside  the  chamber. 


THERMI3TER  METERED 
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Schematic  Diagram  of  the  Resonance  Tube  and  Associated 
Apparatus 
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2.3.2  Impedance  Tube 

The  impedance  tube  is  a  modification  of  the  Bruei  &  Kejer  (B&K) 
standing  wave  apparatus.  A  brass  tube,  100  cm  long  having  inside 
diameter  of  3.81  cm  is  fitted  to  the  speaker  assembly.  At  one  end 
the  threaded  cup  type  of  sealed  fixture  is  rigidly  mounted  to  the  tube; 
see  Figure  4.  This  fixture  accommodates  various  porous  plates  whose 
impedance  is  to  be  determined.  The  tube  and  ancillary  apparatus  are 
shown  in  Figure  5.  The  detailed  description  of  the  equipment  required 
to  drive  the  waves  in  the  tube  and  the  instrumentation  for  measuring  the 
waves  in  the  tube  is  provided  in  Refs.  2-5. 

2.4  Calibration  Procedures  and  Estimates  of  Errors 

The  flow  turning  losses  to  be  determined  are  small  and  the  indirect 
method  of  measuring  them  requires  considerable  care  and  good  instru¬ 
mentation  to  achieve  the  required  accuracy.  References  3  and  4  contain  detailed 
descriptions  of  the  calibration  of  the  gas  supply  system,  extended  here  for 
the  lateral  flow;  calibration  of  microphones  and  sources;  and  estimates  of 
errors  for  both  the  resonance  and  impedance  tubes.  However,  the 
major  source  of  error  in  the  present  work  seems  to  be  associated 
with  measurement  of  the  admittance  functions  of  the  circular  tubes, 
which  required  extrapolation  of  the  admittance  function  data 
for  various  porous  plates.  This  problem  is  discussed  in  the  following 


section. 
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Figure  4.  End  Fixture  for  the  Impedance  Tube 
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ure  5.  Schematic  Diagram  of  the  Impedance 
Apparatus 
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2.  5  Measurement  of  the  Admittance  Functions  for  Cylindrical  Porous 

Tubes 

This  part  of  the  program  became  a  much  larger  problem  than 
we  had  foreseen.  There  is  one  essential  difficulty,  namely  measuring 
the  admittance  or  impedance  of  a  curved  surface.  An  impedance  tube 
is  useful  for  measuring  the  characteristics  of  flat  surfaces  only.  In¬ 
trinsic  constraints  set  by  the  geometry  prevent  direct  measurements  for 
a  carved  surface  by  directing  waves  upon  the  sample  and  measuring 
the  reflections. 

For  application  to  the  apparatus  used  here  we  require  the  response 
of  the  curved  elements  to  pressure  fluctuations  imposed  by  waves 
travelling  parallel  to  the  surface.  The  admittance  function  is  required 
for  conditions  when  mean  flow  is  entering  through  the  surface.  A 
direct  measurement  is  conceivable.  The  sample  may  be  mounted  in  a 
tube  at  one  end  of  which  an  acoustic  driver  excites  travelling  waves. 

The  tube  is  terminated  at  the  other  end  with  an  anechoic  vent  through 
which  both  the  mean  flow  and  the  acoustic  waves  are  exhausted. 
Measurement  of  the  pressure  on  the  upstream  and  downstream  sides 
of  the  porous  sample  should  then  give  the  reduction  of  amplitude  in 
the  travelling  wave  and  therefore  a  direct  determination  of  the  energy 
change  of  the  wave.  The  energy  change  (loss)  can  be  theoretically 
related  to  the  admittance  function  of  the  lateral  surface. 

There  are  several  difficulties,  the  main  one  being  the  very  small 


losses  involved.  The  influence  of  the  surface,  represented  by  its 
admittance  function,  is  not  substantially  different  from  apparent  losses 
caused  by  other  interactions  with  the  non-uniform  flow  field  issuing 
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1 

f  from  the  curved  sample.  That  is,  there  are  necessarily  gradients  of 

the  average  velocity  in  the  axial  direction  along  the  tube.  These  produce 

I 

*  internal  losses  and  scattering  of  the  acoustic  waves.  Furthermore, 

it  is  extraordinarily  difficult  to  construct  an  anechoic  termination  suf¬ 
ficiently  effective  for  our  purposes.  For  reasonable  lengths  of  lateral 
surfaces,  the  changes  in  amplitude  of  the  travelling  waves  are  so  small 
as  to  be  practically  very  difficult  to  measure  with  reproducible  precision 
in  this  kind  of  arrangement.  Because  of  reflections  from  the  termination 
of  the  tube,  waves  travelling  in  both  directions  are  present,  obscuring 
the  data.  And  finally,  the  damning  characteristic  -  the  flow  turning 
loss  is  necessarily  present. 

* 

Consequently,  we  have  no  choice  but  to  measure  the  admittance 
function  for  a  flat  surface  with  waves  incident  in  the  normal  direction 
and  assume  that  the  results  are  applicable  to  a  curve.'  surface  exposed 
to  waves  moving  in  the  parallel  direction.  An  unavoidable  practical 
difficulty,  unforeseen,  then  arose:  we  have  been  unable  after  considerable 
effort  and  time  to  obtain  flat  samples  having  the  same  porosity  -  and 
therefore  flow  characteristics  -  as  the  curved  samples.  It  happens 
to  be  a  consequence  of  the  manufacturing  processes  that  the  cylindrical 
tubes  have  much  lower  porosity  than  any  standard  flat  samples. 

We  have  therefore  made  the  additional  assumption  that  the 
acoustical  properties  of  the  porous  plates  are  defined  uniquely  by 
the  steady-flow  characteristic,  flow  rate  as  a  function  of  pressure 
drop.  The  admittance  functions  for  the  curved  surfaces  (cylindrical 
lateral  flow  elements)  were  determined  by  extrapolation  in  the  following 
way. 

Three  porous  plates,  referred  to  as  grades  F,  G  and  H,  were 

_ 

In  order  to  minimize  costs,  we  used  as  much  as  possible  equipment  from  the 
preceding  work  (Ref.  4);  we  therefore  continued  using  the  circular  resonance 
tube. 
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obtained  from  the  manufacturer,  Pall  Western  Company.  These  were 
first  calibrated,  by  measuring  the  variation  of  flow  rate  with  pressure 
drop  across  a  plate.  The  results,  and  the  results  for  the  two  lateral 
flow  elements  are  shown  in  Figure  6. 

The  admittance  functions  for  the  three  flat  samples  were  then 
measured,  using  an  impedance  tube  as  described  in  Refs.  3-5.  Repeated 
tests  were  made  in  order  to  determine  the  experimental  uncertainty. 

The  data  for  the  real  part  of  the  admittance  function  for  the  three  plates 
are  given  in  Table  1  and  plotted  in  Figures  7-9  for  the  frequency  equal 
to  445  Hz.  The  solid  lines  in  the  figures  are  least  squares  fits  to 
the  data.  Table  2  is  a  summary  of  the  data  for  the  three  flat  samples. 
These  data  were  then  plotted,  as  shown  in  Figure  10,  as  the  real  part 
of  the  admittance  versus  pressure  drop  for  several  values  of  the  average 
Mach  number  of  flow  through  the  plates. 

We  then  make  the  crucial  assumption  that  the  results  may  be 
extrapolated,  as  indicated  in  Figure  10.  The  real  part  of  the  admit¬ 
tance  for  a  cylindrical  element  is  supposed  to  lie  on  the  same  line  of 
constant  Mach  number  (or  flow  speed,  or  flow  rate)  but  at  the  value  of 
pressure  drop  measured  and  shown  in  Figure  6.  Table  3  is  a  summary 
of  values  of  the  real  parts  of  the  admittance  functions  for  the  two 
lateral  flow  elements.  These  results  are  plotted  in  Figures  11  and  12. 
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TABLE  1 


Data  for  the  Admittance  Functions  of  the  Porous  Plate 
Grade  F  (f  =  445  Hz) 


X  104  A^r)  X  102  A^r)  X  102  OX  106 


2. 4984 

0  2.4405  2.4724  240 

2.  4783 

2.  2246 

0.97  2.2379  2.  2230  129 

2. 2064 

1.9224 

1.95  1.9262  1.9103  199 

1. 8822 

2.  2200 

2.92  2.1502  2.  1926  320 

2.  2275 

2.  0549 

3.90  2.0055  2.0237  222 

2.  0106 

2. 2569 

4.87  2.2211  2.2289  204 

2. 2088 
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TABLE  1  (Cont'd.  ) 

Data  for  the  Admittance  Functions  of  the  Porous  Plate 
Grade  G  (f  =  445  Hz) 

X  104  Aj^  X  102  A^  X  102  OX  106 


0. 6719 

0  0. 6079 

0. 6665 

0. 6939 

0.97  0.  6953 

0. 7062 

0. 7060 

1.95  0.6854 

0. 6819 

0. 6784 

2.92  0.6861 

0. 6872 

0. 7083 

3.90  0.6935 

0. 6516 


0. 6488  290 

0.6985  55 

0.6911  106 

0.6839  39 

0. 6845  240 


4.  87 


0. 6295 
0. 6012 
0. 6006 


0. 6104 


135 
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TABLE  1  (Cont'd.  ) 


Data  for  the  Admittance  Functions  of  the  Porous  Plate 
Grade  H  (f  =  445  Hz) 


^  x  104  A^r)  x  102  A^r>  X  102  ax  106 


1.  2504 

0  1.2657  1.2408  252 

1. 2063 

1.  2340 

0.97  1.2344  1.2370  40 

1. 2627 

1. 3958 

1.95  1.3909  1.3860  114 

1.  3702 

1.3283 

2.92  1.3657  1.3538  180 

1.3673 

1. 4445 

3.90  1.4400  1.4380  63 

1.4294 

1.  43  73 

4.87  1.4918  1.4622  225 

1.  4576 


x  10 


ure  8.  Real  Part  of  the  Admittance  Function,  Plate 


150 


Admittance  Function,  Plate 


Summary  of  Results  for  the  Porous  Plates 


Parts  of  the  Admittance  Functions  versus  Pressure  Drop 
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TABLE  3 


* 

* 


Final  Results  for  the  Admittance  Functions 


for 

the  Lateral 

Flow  Element 

CFE  1 

CFE  2 

Flow  Mach  No. 
M  X  104 

A<J>x.02 

a .  x  106 

X 

Flow  Mach  No. 

M  x  104 
t 

X  ‘O2 

at  x  106 

0 

0.  536 

610 

0 

0.  538 

590 

0.97 

0.388 

135 

0.  97 

0.  385 

133 

1.95 

0.  453 

23 

1.95 

0.  449 

23 

2.92 

0.  434 

20 

2.29 

0.  427 

20 

3.  90 

0.  459 

17 

3.  90 

0.  459 

17 

4.87 

0.  363 

11 

4.  87 

0.  360 

11 

Least  squares  fit  to  these  data: 

A(^  =  (.00536  +  .00061)  +  (-  3.277  +  1.  883)M 

tl  ““  ""  *• 

a£>  =  (.00538  +  .00059)  +  (-  3.143  +  1.  825)  Mfc 


I 
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2.  6  Measurements  of  the  Flow  Turning  Loss 

With  the  results  of  Refs.  3-5  and  the  preceding  section,  we  have 

the  values  of  a  ,  a,  and  a  in  eq.  (2.32),  all  as  functions  of  the 
p  t  v 

appropriate  Mach  number.  Before  can  be  measured  it  is  necessary 

to  measure  the  contribution  from  other  losses,  denoted  by  a^.  The 
series  of  tests  carried  out  to  achieve  this  result  also  serve  as  an 
important  check  of  the  experimental  procedures  by  providing  a  means 
of  assessing  the  self-consistency  of  the  data.  It  is  during  this  portion 
of  the  experimental  work  that  many  faults  with  the  apparatus  have  been 
found  and  corrected  (see  also  Refs.  3  and  4).  We  discuss  here  only 
the  final  series  of  tests. 

Table  4  contains  data  taken  for  three  series  of  tests  performed 
at  each  Mach  number  M^  for  the  average  flow  through  the  lateral 
boundary.  The  average  values  and  standard  deviations  of  both  the 
net  attenuation  constant  and  resonance  frequency  are  also  given. 

Figure  13  shows  the  variation  of  the  average  values  of  the  net  attenua¬ 
tion  constant  with  respect  to  the  mean  cross-flow  Mach  number  M^ 
at  the  resonance  frequency  near  445  Hz.  The  solid  line  in  the  figure 
represents  the  least  squares  fit,  assuming  that  Qnet  varies  linearly 
with  M^. 

Two  series  of  tests,  as  shown  in  Table  4,  were  carried  out 
with  no  flow:  one  series  with  tne  vent  open  and  one  with  the  vent 
closed.  The  difference  in  the  attenuation  constants  should  represent 
the  effects  of  radiation  losses.  Here,  the  difference  is  -  0.01  + 

0.02  sec  *,  representing  a  small  average  gain,  although  it  could 
(as  it  should  be)  be  a  loss  within  the  experimental  uncertainty. 


L 


J 
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TABLE  4 

Data  for  the  Attenuation  Coefficient  with  Lateral  Flow 
(1/2"  vent,  f  =  445  Hz.) 


Mpx  104 

Mt  x  104 

anet 
(sec- *) 

f 

Hz 

a 

net 

(sec-  1 ) 

a 

a 

(sec- * ) 

f 

Hz 

°f 

Hz 

0 

0 

63.99 
64.27 
64.3  6 

446. 00 
446.  60 
446. 50 

64.21 

0.  16 

446.  40 

0.28 

4.  43 

0 

64.22 

64.  24 
64.  30 

445. 90 
446. 00 
446. 00 

64.  25 

0.03 

445.  97 

0.  05 

4.  43 

0.97 

64.  69 
65.90 

66.  03 

445.  40 
445. 50 

446.  50 

65.  87 

0.  14 

445.  80 

0.  50 

4.  43 

1.95 

66.  21 

66.  69 

66.  61 

446.  20 

447.  90 
447. 00 

66.  50 

0.  21 

447. 03 

0.69 

4.  43 

2.92 

66.  65 

66.  60 
66.35 

445.00 
445. 50 
445. 50 

66.  53 

0.  13 

445.33 

0.24 

4.43 

3.  90 

66.  79 

67.  08 
66.95 

446. 00 
446. 50 
446. 50 

66.  94 

0.  12 

446. 33 

0.24 

64.04  445.00 

0  0  64.25  445.50  64.22  0.13  445.33  0.24 

(vent  closed)  64.35  445.50 
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Because  of  the  experimental  difficulties  discussed  at  greater  length  in 
Refs.  3  and  4,  the  numerical  value  of  the  measured  radiation  loss  is 
necessarily  uncertain. 

The  consistency  of  the  data  and  losses  through  the  lateral  boundary 
have  been  determined  with  no  flow  and  the  vent  closed,  using  the  equation 

a  =  a  +  a  +  a  (2.3  5) 

p  t  d 

The  value  of  a  was  measured,  and  with  a  and  a,  given  by  pre- 

p  t 

vious  results,  we  have  found 


a,  =  14.90  +  2.75  sec 
a  — 


-1 


(2.  36) 


The  theoretical  value  for  viscous  and  heat  transfer  losses  calculated 
according  to  classical  formulas  for  laminar  flow  is  11.20  sec  The 

difference,  3.70  +  2.75  sec" 1,  represents  other  contributions,  probably 
located  mainly  at  the  driver  pistons,  with  a  small  radiation  loss  through 
the  lateral  boundary.  The  relatively  small  size  of  these  other  con¬ 
tributions  is  a  result  of  the  care  taken  to  perfect  the  apparatus. 

According  to  the  analysis  described  in  Section  2.  1,  the  decay 
constant  associated  with  the  flow  turning  losses  is  expected  to  vary 
linearly  with  the  Mach  number  of  the  average  flow  entering  the  boundary. 
For  the  configuration  shown  in  Figure  2,  the  formula  (2.  28)  gives 


‘ft  -  2fMt  Xs 
c 


(2.37) 


where  A.  is  a  geometrical  factor,  evaluated  in  Refs.  13  and  17,  which 

Xj 

is  equal  here  to  0.97.  We  therefore  assume,  in  correlating  our  data, 
that  a^t  ~  Mj.  and  we  shall  compare  the  slope  da^/dM^  computed  for 
the  data  with  that  given  by  (2.3  7). 


We  assume  that  the  Mach  numbers  and  frequency  are  given  with 

no  error,  so  that  eq.  (2.32)  can  be  written  to  show  explicitly  the  average 

values  and  standard  deviations,  after  (2.33)  and  (2.34)  have  been  used: 

af.  =  +  Aa,,  =  (a  +  Aa)  -  [4f(A^  +  AA^  +  M  ) 

ft  ft  —  ft  —  '  d  —  p  p 


+  4T(A<r)  t  4  A<r)  +  M  ,  JiS  O 

c 

+  K  ±  >  +  (“ri  i 


i 


(2.38) 


— .  /  j  ( r )  — { r )  ( r ) 

The  values  for  A  ,  A  A  ,  A.  ,  AA  are  averages  for  the  components 
p  p  t  t 

in  the  two  sides  of  the  resonance  tube.  All  of  the  quantities  in  the  square 
brackets  on  the  right  hand  side  of  (2.  38)  are  known  from  the  previous 
experimental  work  described  above.  These  are  summarized  in  Tables 
5-7  for  direct  use  in  (2.38). 

Equation  (2.38)  may  now  be  used  with  the  data  summarized  in 
Tables  4-7  to  infer  the  values  of  a^.  These  are  given  in  Table  8 
and  plotted  in  Figure  14.  The  solid  line  in  Figure  14  represents  a 
least  squares  fit  to  the  data  using  the  equation 


—  daft 

aft  =  (a,. ) _  +  M  {-2-  ) 

Mt=0  dMt 


(2. 39) 


The  least  squares  fit  was  done  without  imposing  the  physical  requirement 
that  the  linear  fit  should  pass  through  origin,  using  a  library  routine  in 
its  primitive  form.  Note  that  the  experimental  uncertainty  in  at 
M^=0  encompasses  the  average  of  the  measured  values  is  slightly 

negative. 

The  derivative  dav/dMt  of  equation  (2.38)  can  be  used  with  the 
date  to  provide  an  independent  result: 
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TABLE  5  Values  of  the  Admittance  Function  for  the  Faces  of  the 
Driver  Pistons,  f  =  442  Hz.  (Ref.  4) 


M 

A<r> 

0 

P 

P 

0  -4 

.01544 

. 00009 

4.43  x  10 

. 01550 

.  00019 

TABLE  6  Values  of  the 

Attenuation 

Coefficient  for  the 

Vent,  f  =  442 

Hz.  (Ref. 

4) 

M  x  104 

a 

Aq 

V 

V 

V 

0 

0.06 

0.  70 

4.  43 

-0.  59 

0.  85 

8.  86 

-0.90 

0.  78 

13.  29 

-1.  15 

0.  72 

17.  72 

-2.  39 

1. 00 

22.  15 

-3.  15 

0.  73 

TABLE  7  Values  of  the 

Admittance 

Function 

for  the  L 

Elements,  f  = 

445  Hz  (Average  of 

values  in 

Mt  X  104 

A(r) 

At 

a 

0 

0. 5320 

.  0006 

0.97 

0.  3862 

.  0013 

1.95 

0. 4512 

.  0002 

2.92 

0.  4307 

.  0002 

3.  90 

0.  4576 

.  0002 

< 


ure  M.  Measured  Attenuation  Coefficient  for  Flow  Turning  Losses 
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daft  daft 

+  A  (-— 

dMfc  dM 


)  =  —7-  +  A  (-%■) 
dM  dM 

dA;r)  dA[r) 

.  4f  {— A_  +  £  )  +  1} 


dM, 


dM, 


da  da  dM 

-  K-dr>  ±  * 

dM  dM  dM. 

v  t 


(2.40) 


The  value  using  the  data  discussed  above  is 

d“ft  daft  - 1 

+  A  (-~)  =  23000  +  11000  sec  (2.41) 

dMt  dMt 

This  is  five  times  as  large  as  the  theoretical  value  4000  sec  * 
calculated  using  eq.  (2.3  7). 

There  is  little  doubt  that  there  is  indeed  a  loss  of  acoustic  energy 
associated  with  the  flow  in  at  the  lateral  boundary,  and  that  it  increases 
with  flow  rate,  or  Mach  number.  Unfortunately,  owing  to  the  large 
experimental  uncertainty  associated  with  these  results,  it  is  not  possible 
to  state  with  any  degree  of  confidence  how  closely  the  one-dimensional 
approximation  represents  the  actual  losses. 

We  emphasize  in  conclusion  that  there  are  two  major  sources  of 
uncertainty:  the  necessity  to  obtain  the  admittance  function  for  the 
curved  lateral  flow  elements  by  extrapolation  from  the  data  taken  for 
flat  elements;  and  the  problem  of  obtaining  the  desired  quantities  as 
small  differences  of  large  numbers.  The  only  way  to  overcome  the 
second  difficulty  is  to  use  a  larger  scale  apparatus,  requiring  a  larger 
gas  supply  not  available  to  us.  If  a  sufficiently  large,  high  pressure 
supply  is  available,  all  porous  elements  should  be  operated  choked  to 
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prevent  spurious  coupling  between  the  acoustic  field  in  the  chamber  and 
the  system  upstream  of  the  porous  elements.  We  believe  that  those 
losses  have  been  minimized  in  our  tests,  although  we  have  been  able 
to  provide  choked  flow  only  through  orifices  further  upstream. 
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III.  MEASUREMENTS  OF  THE  ADMITTANCE  FUNCTION  OF  A 

CHOKED  NOZZLE 

The  objective  of  this  portion  of  the  program  was  to  develop  a  technique 
for  measuring  the  admittance  function  of  a  choked  nozzle,  including  the  effects 
of  mean  flow  and  nozzle  submergence.  The  approach  taken  involved  adaptation 
of  a  continuously  operating  supersonic  tunnel.  This  approach  was  expected  to 
have  several  advantages  over  use  of  the  blow-down  facility  used  in  the  most 
recent  previous  work  (Refs.  9-12):  the  capitol  investment  in  measuring 
equipment  and  instrumentation  is  considerably  less;  it  is  easier  to  study  the 
flow  field,  for  example,  to  determine  whether  the  field  is  really  one-dimensional 
as  required  by  the  usual  theory;  reduction  of  data  is  often  simpler;  and  tests 
to  obtain  comparable  information  should  be  less  expensive  to  conduct. 

During  the  early  part  of  the  work,  we  carried  out  tests  to  deter¬ 
mine  the  severity  of  contamination  of  the  test  section  by  noise  and  oil 
from  the  compressors  running  the  tunnel.  The  amount  of  oil  present 
was  a  serious  problem,  prompting  the  decision  to  operate  the  tunnel 
as  an  open  circuit.  A  diagram  of  the  entire  facility  is  given  in  Figure 
15,  and  the  open-circuited  test  section  is  shown  in  Figure  16.  Air  is 
drawn  in  from  the  room  and  flows  through  the  test  section  which  ter¬ 
minates  in  the  choked  nozzle  to  be  tested. 

Acoustic  waves  are  generated  at  the  entrance  to  the  test  section 
by  an  acoustic  driver,  an  Atlas  sound  driver.  Model  5D-370.  The  test 
section  then  serves  as  an  impedance  tube  and  the  primary  measurements 
are  taken  to  determine  the  distributions  of  amplitude  and  phase  for  the 
standing  acoustic  field  within  the  test  section.  In  principle,  the  tech¬ 
nique  is  sound  and  is  essentially  that  used  previously  in  Refs.  9-12. 

In  practice  there  are  substantial  difficulties  which  we  had  not  foreseen. 

It  appears  that  they  can  be  overcome  for  a  broad  range  of  conditions  and  the 


FLARED 
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;ure  16.  Sketch  of  the  Test  Section  in  the  Supersonic  Wind  Tunnel 
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method  should  be  regarded  as  a  promising  technique  for  obtaining  the 
admittance  functions  for  choked  nozzles.  There  are  some  important 
improvements  which  are  evident  from  the  present  work;  these  will  be 
noted  in  the  following  sections  and  summarized  in  the  concluding  remarks. 

3.  1  Analysis  of  the  Impedance  Tube 

The  traditional  applications  of  the  impedance  tube  have  two 
important  characteristics;  there  is  no  average  flow  through  the  tube; 
and  the  loss  of  acoustic  energy  at  the  test  sample  are  generally  quite 
large,  producing  relatively  large  changes  in  the  pressure  amplitude 
in  the  acoustic  field.  Both  of  those  characteristics  are  significantly 
changed  in  the  present  application. 

Flow  through  the  tube  produces  a  boundary  layer  which  in  fact 
becomes  quite  thick  for  the  lengths  of  tubes  required- 
This  must  have  some  influence  on  the  acoustic  field,  but  we  have  not 
examined  this  problem  in  this  work.  We  have  assumed  throughout  that 
the  acoustic  field  is  one-dimensional,  with  phase  fronts  flat  and  normal 
to  the  axis  of  the  tube.  Some  limited  tests,  described  later  in  Section 
3.3,  have  shown  that  the  assumption  is  quite  good,  for  the  apparatus 
used  in  these  experiments,  if  the  frequency  is  less  than  1000  Hz.  and 
the  Mach  number  of  the  average  flow  is  less  than  0.  20. 

The  procedures  eventually  used  for  reducing  the  data  are  described 
in  Section  3.4.  Here  we  describe  the  simplest  analysis  of  the  im¬ 
pedance  tube  with  flow.  As  sketched  in  Figure  17,  the  stationary 
acoustic  field  may  be  represented  as  the  superposition  of  two  travelling 
waves.  Let  the  nozzle  be  located  at  z  =  0,  with  the  test  section  stretching  over 


z  >  0.  With  a  uniform  flow,  in  the  negative  z-direction,  the  governing  equations 
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are  (2.11)  and  (2.12)  without  the  sources: 


—  9u'  + - du'  9p' 

"  sr  pu  ST  +  *  0 


(3.  1) 


^El  +  U  3£l  +  y  P  9^  =  o 

at  u  az  rp  az  u 


(3.2) 


The  field  in  the  tube  consists  of  a  wave  travelling  to  the  right,  incident 
upon  the  nozzle,  and  the  reflected  wave  travelling  to  the  left;  in  steady 
state  the  superposition  of  the  two  waves  forms  a  stationary  wave  pattern. 
Appropriate  solutions  to  (3.1)  and  (3.  2)  are 


ik+z  -ik  z  -iuot 

p'  =  [P+  e  +  P_  e  ]  e 


(3.3) 


ik  z  -ik  z  -iuut 

u'  =  [U+  e  +  U  e  ]  e 


(3.4) 


Waves  to  the  right  (left)  are  denoted  by  +  (-).  Substitution  of  (3.3) 
and  (3.4)  into  (3.1)  and  (3.2)  leads  to  four  equations  for  the  coefficients 
P+J  P  ,  U  ,  U  .  The  condition  for  non-trivial  solutions  gives  the 
formulas 


1  -  M 


1+  M 


(3.  5) 


where  here  k  =  u)/a  is  the  real  wavenumber  because  no  losses  have 
been  accounted  for.  The  equations  may  then  be  used  to  relate  the 
coefficients; 


U  =  — p  ; 

+ - +  ’ 

P  a 


U  =  -J—  P 
P  a 


(3.6) 


Then  (3.3)  and  (3.4)  may  be  written 


p'  =  [P+  e  +  P_  e  ]  e 


■iKz  -i(uut  -  MKz) 


(3.  7) 


-67- 


,  iKz  -iKz  -i(iut-MKz) 

u'  =  —  [P+  e  -  P_  e  ]  e 
P  a 


(3.  8) 


1  -  M2 


(3.9) 


The  admittance  function  of  the  nozzle  is  defined  as 


- u1 

A  =  p  a 


(3.  10) 


Substitution  of  (3.7)  and  (3.8)  into  this  formula  leads  to  the  relation 


between  P  and  P+: 


p+  =  5P.  s  6  =  “  l6le 


(3.  11) 


The  reflection  coefficient,  6,  has  amplitude  and  phase  given  by  the 


formulas 


1  -  A 


-  2AV 


(1+aM)2  7  (A(i»,2  1  =  ,1+A<r>>*  +  (A*1')2 


(3.  12) 


The  pressure  field  may  be  written  in  the  form 

iKz  -iKz  -i(out  -  MKz) 

P'  =  P_  [e  +  6e  ]  e 


=  Pe 


-i(u)t  -  MKz  -  \jf 


where,  with  P  =1, 

P  =  1  +  1  6 1  2  +  Z  J  6 1  cos  (ZKz  +  cp) 


(3.  13) 


(3.  14) 


ten  *  =  -  sinKz  +  l6lsin(Kz  +  cp) 
™  cosKz  +  6  cos  (Kz  +  cp) 


(3.  15) 


Equation  (3.  13)  shows  explicitly  how  the  pressure  field  may  be 
written  in  terms  of  its  amplitude  and  phase.  A  similar  form  may  be 


constructed  with  viscous  losses  included,  but  there  is  no  reason  here 
to  do  so.  The  purpose  of  the  preceding  is  to  display  as  simply  as 
possible  the  structure  of  the  acoustic  field,  particularly  the  way  in 
which  the  amplitude  and  phase  depend  on  the  admittance  function. 

The  experimental  procedure  is  based  on  the  idea  that  data  for  the 
amplitude  and  phase  can  be  interpreted  with  an  analytical  result,  such 
as  (3.  13)  to  give  the  real  and  imaginary  parts  of  the  admittance  function. 
The  details  will  be  pursued  further  in  Section  3.  4  where  the  method  of 
data  reduction  is  discussed. 

3.2  Description  of  GALCIT  Supersonic  Wind  Tunnel  and  Experimental 


The  schematic  of  the  GALCIT  (Guggenheim  Aeronautical  Laboratory, 
California  Institute  of  Technology)  supersonic  wind  tunnel  is  given  in 
Figure  14.  The  tunnel  is  designed  to  operate  continuously.  The  com¬ 
pressors  (maximum  compression  ratio  of  3.  2  each)  were  run  in  series 
with  a  200  HP  electric  motor.  In  the  present  work,  one  of  the  test 
sections  was  disconnected  from  the  main  flow  line  and  the  other  was 
run  open  circuit.  The  open  circuit  arrangement  makes  it  possible  to 
use  the  room  air  which,  after  expanding  through  the  nozzle,  is  compressed 
by  the  compressors  and  then  exhausted  back  into  the  atmosphere.  This 
passage  of  air  helps  reduce  the  flow  noise  otherwise  transmitted  to  the 
test  section  from  the  running  of  the  compressors.  Also,  it  prevents 
the  deposition  of  the  compressor  lubricating  oil  to  the  surface  of  the 
microphones. 

The  stagnation  pressure  in  the  test  section  is  therefore  maintained 
to  be  atmospheric.  The  compressed  air  from  the  compressor  passes 
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through  an  aftercooler,  one  for  each  compressor.  The  heat  from  this 
compressed  air  is  exhausted  by  cooling  water  which  in  turn  is  cooled 
by  the  cooling  tower  on  the  roof  of  the  laboratory.  To  maintain  the 
compression  ratio  within  the  desired  limits,  the  valve  settings  for  the 
various  valves  were  adjusted  for  each  flow  Mach  number.  With  appropriate 
cooling,  test  run  times  of  approximately  2  hours  were  easily  obtained. 

The  impedance  tube  used  to  measure  the  admittance  of  the  nozzles 
is  sketched  in  Figure  16.  The  test  section  is  an  aluminum  tube  150  cm 
long  having  an  internal  diameter  of  5  cm  and  wall  thickness  of  0.64  cm. 

At  one  end,  a  radiating  horn  shell  is  connected  which  serves  as  the 
flared  intake  for  incoming  air,  while  at  the  other  end  a  nozzle  is 
mounted.  The  flow  through  the  nozzle  is  carried  to  the  suction  end 
of  the  compressors  for  the  supersonic  wind  tunnel.  A  1.25  cm  B&K 
microphone.  Model  4133,  is  mounted  on  the  tube  wall  at  a  distance  of 
10  cm  from  the  entrance  plane  of  the  nozzle.  Another  microphone  is 
mounted  on  elastic  supports  within  the  microphone  trolley.  The  micro¬ 
phone  probe  is  275  cm  long,  open  at  one  end  and  connected  to  the  micro¬ 
phone  trolley  at  the  other  end. 

The  probe  is  supported  at  four  different  locations  within  the  tube 
using  teflon  tipped  spider  supports.  The  trolley  slides  on  a  graduated 
scale  fitted  with  its  own  bernier.  The  major  advantage  of  this 
arrangement  is  the  accuracy  with  which  the  acoustic  minima  can  be 
located.  The  trolley  is  filled  with  cotton  balls  to  aid  irsulation  against 
both  external  air-born  noise  and  structural  vibrations.  Near  the  flare 
intake,  a  100  Watt  Atlas  sound  acoustic  driver  Model  SD-370  is 
mounted  on  the  tube  wall.  This  driver  imposes  acoustic  oscillations 
on  the  main  flow  of  air. 

ti 
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Table  9  Summary  of  the  Characteristics  of  Nozzles  Used  for 
Admittance  Measurements 


Mach  Number 

Secondary 
Flow,  $  of 

Throat 

Diameter, 

d* 

Submersion 

of  the  Main  Flow 

Main  Flow 

(cm) 

(cm) 

0.  04 

-- 

1.  25  > 

0.  10 

-- 

2.  11 

0.20 

-- 

2.95 

>- 

small 

conical 

0.  30 

3.66 

nozzles  in 

0.35 

3.91  v 

plates 

0.  04 

0.  04 

0.  04 


0.  48 

0 

1.25 

4.  29 

0.48 

2.  5 

1. 25 

4.29 

0.48 

5.0 

1.  25 

4.  29 
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Table  9  is  a  list  of  the  different  nozzles  constructed  for  the 
experiments.  One-dimensional  compressible  flow  formulas  were  used 
to  determine  the  throat  diameters  of  the  nozzles.  For  simplicity,  both 
the  convergent  and  divergent  sections  of  the  nozzles  were  conical. 

Figure  18.  The  nozzles  are  machined  of  aluminum.  Only  one  of  these 
nozzles  was  constructed  to  accommodate  the  secondary  flow  injection; 
see  Figure  18.  A  grade  "H"  porous  sintered  steel  plate  (Pall  Western 
Co.)  0.32  cm  thick  was  glued  to  the  front  end  of  the  nozzle  manifold. 

This  nozzle  was  fabricated  with  external  threads.  This  arrangement  with 
additional  aluminum  spacers  permits  variable  submersion  of  the  nozzle 
in  the  main  flow. 

In  addition  to  the  main  air  supply,  the  apparatus  and  the  ancillary 
equipments  shown  in  Figure  19  can  be  divided  into  three  subsystems: 

(i)  the  equipment  required  to  drive  the  acoustic  waves  in 
the  tube 

(ii)  instrumentation  for  measuring  the  waves  in  the  tube 

(iii)  secondary  air  flow  supply  and  instrumentation  necessary 
to  measure  and  monitor  the  flow. 

A  Wavetek  frequency  synthesizer.  Model  171,  together  with  a 
McIntosh  amplifier.  Model  MI-75,  was  used  to  excite  the  driver.  The 
Mura  power  meter.  Model  APMG-18,  was  connected  in  series  to  monitor 
the  electrical  power  input  to  the  acoustic  driver.  A  Monsanto  counter 
timer.  Model  1 0 1  B,  was  used  to  monitor  the  output  frequency  of  the 
synthesizer.  Pressure  signals  from  the  end  of  the  probe  tube  at  the 
wall  were  detected  by  1.25  cm  B&K  microphones.  Model  4133,  as 
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Figure  18.  Drawings  of  the  Nozzle  Tested 


FLOW  IN 


Figure  19.  Instrumentation  for  Measurements  of  the  Nozzle 
Admittance  Function 
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described  earlier.  A  track  scale  was  used  to  measure  the  distance  of 
the  location  of  the  front  end  of  the  probe  tube  from  the  nozzle  face. 

For  a  large  majority  of  the  experiments,  both  signals  were  filtered 
through  a  two  channel  unity  gain  tracking  filter,  Spectral  Dynamics  Model 
SD-122,  with  relative  phase  of  0°  and  1  Hz.  band  pass  filters  driven  by  a 
carrier  generator.  Spectral  Dynamics  Model  SD-120.  The  filtered  output 
of  the  microphone  in  the  trolley  was  amplified  using  a  Hewlett-Packard 
amplifier.  Model  465  A  and  was  then  used  to  measure  the  relative  phase  with 
respect  to  the  wall  microphones  filtered  signal  using  a  Wavetek  phase 
meter.  Model  740.  For  some  of  the  experiments  the  desired  signal  was 
extracted  from  the  noise  by  using  a  Princeton  Applied  Research  Lock-in 
Amplifier,  Model  124  A  with  preamplifier.  Model  117.  The  lock-in 
amplifier  was  externally  tuned  to  a  filtered  signal  from  the  wall  microphone 
using  a  Princeton  Applied  Research  selective  amplifier.  Model  189.  For 
some  experiments,  the  signal  from  the  frequency  synthesizer  was  used 
to  tune  the  lock-in  amplifier.  The  former  technique  of  conditioning  the 
signals  using  the  tracking  filter  enabled  us  to  measure  the  actual  signal 
to  noise  ratio  near  the  acoustic  pressure  minima.  That  information  cannot 
successfully  be  obtained  using  the  lock-in  amplifier.  The  output  of  the 
lock-in  amplifier  and  all  other  AC  and  DC  outputs  used  for  monitoring 
various  parameters  are  measured  with  a  5-1/  "igit  Data  Precision 
Voltmeter,  Model  3  500. 

The  calibration  of  the  secondary  air  flow  supply  was  done  as 
described  in  Ref.  3.  The  pressure  drop  across  the  porous  plate  was 
monitored  with  a  Statham  differential  pressure  gauge  (0-50  psi).  A 
Dynamic  bridge  supply.  Model  6343  A  was  used  with  the  pressure  gauge. 


-75- 


The  porous  plate  was  chosen  to  provide  very  high  acoustic  impedance, 
which  implies  a  relatively  large  pressure  drop  for  a  given  secondary 
flow  rate. 


3.  3  Preliminary  Tests 

A  large  number  of  preliminary  tests  were  conducted  during  the 
early  phase  of  this  program  to  characterize 

(a)  the  flow  noise 

(b)  acoustic  pressure  symmetry  in  the  tube,  and 

(c)  the  influence  of  the  support  mount  of  the  long  probe  on 
the  acoustic  pressure  measurements. 

Results  from  these  experiments  were  used  to  guide  the  selection  of 
appropriate  parameters  such  as  the  flow  Mach  number  and  the  strength 
of  the  acoustic  pressure  signal  during  the  tests. 

An  earlier  test  was  done  with  a  thick  plate  orifice;  the  Mach  number 
of  the  mean  flow  was  approximately  0.  20.  The  acoustic  signal  was 
measured  at  a  fixed  location  from  the  orifice.  It  was  found  that  the 
flow  noise  was  102  db  +  2  db  (reference  pressure  2  X  10  N/m  )  over 
a  frequency  range  of  300  to  3000  Hz.  It  was  subsequently  determined 
that  the  flow  noise  is  independent  of  the  location  of  the  probe.  Numerous 
tests  were  conducted  at  flow  Mach  numbers  of  0.  10  or  less.  Those 
tests  verified  that  the  signal  to  noise  ratio  at  the  acoustic  pressure 
minima  was  at  least  2.  8  db.  This  figure  was  measured  and  checked 
using  a  finely  tuned  tracking  filter  as  described  earlier. 

Non-uniformities  in  the  acoustic  pressure  over  the  cross-section 


of  the  tube  may  arise  because  of  the  presence  of  a  non-uniform  flow 
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distribution  or  because  of  the  location  of  the  sound  source  which  is  in 
the  tube  wall.  A  check  of  the  symmetry  of  the  acoustic  pressure  was 
considered  essential  for  subsequent  treatment  of  these  data  using  the 
one-dimensional  theory  of  an  impedance  tube.  Some  experiments  were 
carried  out  to  determine  whether  or  not  the  phase  fronts  are  flat.  Two 
microphones  were  placed  exactly  opposite  each  other  ten  centimeters 
from  the  thick  plate  orifice  (flow  Mach  number  of  approximately  0.  20). 
There  appears  to  be  negligible  asymmetry  for  frequencies  below  750  Hz.  , 
but  at  2000  Hz.  a  difference  of  about  20$  was  found  between  the  ampli¬ 
tudes  on  opposite  sides  of  the  tube. 

For  the  2-inch  diameter  tube  used  here,  the  frequency  of  the 
lowest  tangential  mode  is  roughly  3  500  Hz.,  so  transverse  modes  are 
not  excited  directly  by  the  acoustic  driver.  The  asymmetry  in  the  field 
is  due  to  the  location  of  the  driver  on  one  side  of  the  tube.  We  were 
unsuccessful  in  some  tests  using  two  opposed  drivers.  The  problem  of 
obtaining  a  proper  match  of  both  amplitude  and  phase  of  the  two  drivers 
is  a  very  difficult  one  and  the  fields  we  obtained  in  this  way  were  not 
as  uniform  as  those  excited  by  a  single  driver.  A  single  driver  mounted 
on  the  axis  at  the  inlet  is  the  best  arrangement,  but  that  requires  special 
hardware  when  the  pressure  probe  is  also  located  on  the  axis,  as  we  chose 
here. 

As  discussed  earlier,  a  long  probe  was  used  to  measure  the  acoustic 
pressure  distribution  along  the  length  of  the  tube.  A  large  number  of 
measurements  are  easily  obtained,  distributed  along  the  axis  in  any 
desired  manner.  This  improves  the  accuracy  of  the  results  obtained 
from  the  numerical  scheme  used  to  reduce  the  data.  The  long  probe 
was  supported  at  three  to  five  different  locations  using  sprocket  type 


-77- 


supports.  These  supports  move  together  with  the  probe  tube  sliding 
along  the  tube  walls  during  measurements.  Their  locations,  shape, 
size  and  the  presence  of  the  probe  is  bound  to  affect  the  acoustic 
pressure  signal  at  the  measuring  positions  downstream.  The  critical 
measure  of  this  influence  was  characterized  using  the  solid  end  plate 
and  the  condition  of  no  flow  in  the  tube.  Two  microphones,  one  at  a 
distance  10  cm  and  the  other  at  a  distance  of  20  cm  from  the  solid 
end  plate,  were  mounted  in  the  wall  of  the  tube.  The  test  frequency 
was  selected  to  be  425  Hz.  so  that  the  first  acoustic  pressure  minimum 
was  at  the  second  probe,  20  cm  from  the  end  plate.  For  a  fixed  source 
excitation,  the  amplitude  ratio  and  the  phase  difference  between  these 
two  microphones  were  recorded  as  a  function  of  the  location  of  the  probe. 
It  was  found  that  when  the  probe  is  at  a  distance  of  more  than  25  cm 
from  the  end,  the  amplitude  ratio  and  the  phase  difference,  approximately 
45°,  are  not  affected  by  its  presence.  If  the  support  is  close  to  the  end 
of  the  probe,  both  the  amplitude  ratio  and  phase  are  strongly  affected 
when  the  probe  is  near  the  end  wall. 

This  can  be  attributed  to  the  presence  of  the  sprocket  support.  In 
all  the  tests  reported  here,  the  support  was  therefore  located  at  least 
15  cm.  from  the  end  of  the  probe.  This  configuration  minimizes  the 
influence  of  the  support  and  we  have  made  no  corrections  to  the  data 
to  account  for  the  presence  of  the  probe  and  its  support. 

3.  4  Procedures  for  Data  Reduction 

When  this  program  began,  we  anticipated  that  there  should  be  no 
difficulty  using  the  well-known  procedure  for  determining  the  admittance 
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function  from  measurements  of  the  pressure  amplitude.  The  procedure 
is  straightforward  and  had  apparently  been  used  successfully  in  the  work 
reported  in  Refs.  9-12.  However,  we  encountered  severe  difficulties 
which  consumed  considerable  time  before  we  realized  their  source. 

In  summary,  the  main  conclusion  is  that  measurements  of  the 
pressure  amplitude  only  are  adequate  to  determine  accurately  the  real 
and  imaginary  parts  of  the  admittance  function  providing  the  real  part 
is  relatively  large.  That  is,  the  losses  of  acoustic  energy  incurred  at 
the  end  of  the  tube  must  be  sufficiently  large  -  therefore  making  the  real 
part  of  the  admittance  large  -  that  the  structure  of  the  acoustic  field  is 
substantially  perturbed  from  the  ideal  case.  The  losses  at  the  end 
affect  most  significantly  the  structure  of  the  field  in  the  vicinity  of  the 
pressure  minima,  where  the  experimental  error  is  greatest. 

If,  as  happens  to  be  the  case  for  the  influence  of  the  choked  nozzles 
used  in  the  present  work,  the  losses  are  not  large  enough,  the  traditional 
procedure  leads  to  large  errors  and  in  many  cases  results  having  the 
wrong  sign  (i.  e.  a  gain  of  energy  may  be  inferred).  Under  these  cir¬ 
cumstances,  it  is  absolutely  essential  that  both  the  amplitude  and  the 
phase  of  the  pressure  field  be  measured  and  that  both  pieces  of  informa¬ 
tion  be  used  in  the  data  reduction. 

The  importance  of  phase  measurements  was  earlier  recognized  by 
Baum  (Ref.  18)  and  was  brought  to  our  attention  while  the  present  program 
was  in  progress.  We  had  encountered  difficulties  while  reducing  data 

18.  Baum,  J.  D.  ,  Daniel,  B.  R.  and  Zinn,  B.  T.  "Determination  of  Solid 
Propellant  Admittances  by  the  Impedance  Tube  Method",  AIAA  18th 
Aerospace  Sciences  Meeting  (Jan.  1980)  AlAA  Paper  80-0281. 
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take  n  primarily  to  check  the  apparatus,  experimental  technique  and  the 
procedures  for  data  reduction.  We  then  spent  considerable  effort  ex¬ 
amining  the  sensitivity  of  the  traditional  method  using  "data"  calculated 
for  idealized  cases.  It  is  quite  clear  now  that  a  technique  based  on  data 
for  both  amplitude  and  phase  should  always  be  used.  We  have  adopted 
without  change  the  numerical  procedures  suggested  in  Ref.  18  and  have 
found  them  to  be  very  accurate  and  easy  to  use. 


3.  4.  1  Application  of  Analysis  of  the  Impedance  Tube 

Before  summarizing  the  techniques  finally  used,  it  is  worthwhile 
demonstrating  what  goes  wrong  with  the  traditional  method,  the  use  of 
which  can  no  longer  really  be  justified. 

With  a  uniform  flow  field,  and  distributed  viscous  losses  along 
the  tube  accounted  for,  the  leftward  and  rightward  running  waves  sketched 
in  Figure  16  and  introduced  in  Section  3.  1  can  be  written 

-tCz  i(k  z  -  a't)  1  -Kz  i(k  z  -  uut) 

p^  =  P+e  e  u'  =  -  P  e  e 

P  a 


P' 


<z  -(k  z  +  wt) 
Pee”  u' 


.  Kz  -i(k  z+uut) 

Pee 
P  a 


(3.  1 6)a, b 


The  b'  dary  condition  at  z  =  0  is  satisfied  in  the  form 
P+  =  -  P_e  ;  \|r  =  ttq0  +  iir(30 


and  the  stationary  pressure  field  is  the  superposition  of  p^  and  p'  : 
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Now  define 

■tr(a  +  if3)  =  -  Kz  +  iKz  +  ty  =  (ttQq-Kz)  +  ifirPg+Kz) 

and  the  pressure  field  is 
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The  amplitude  and  phase  of  the  spatial  distributions  are 

|  p  |  =  2P  e  0  [cosh2  ^(a,.  -  — )  -  cos2  tt(P  +  — — ■  )] 

U  11  U  Ml-M4) 


1/2 
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0  =  +  Kz  +  tan  *  [cothir(an  -  —  )tan  Tr(Pn  + — — )]  (3.21) 

0  0  TT  0  X(1_m2} 

where  the  wavelength  is  \  =  2ir/k. 

The  formula  (3.21)  is  the  basis  for  the  familiar  method  of  treating 
data  taken  with  an  impedance  tube,  as  for  example,  explained  in  Ref.  19 
There  are  four  parameters  in  eq.  (3.21):  tig,  Pq,  K,  and  the  amplitude 
P  .  The  real  and  imaginary  parts  of  the  admittance  function  are  found 


19.  Morse,  P.  Vibration  and  Sound.  McGraw-Hill  Book  Co.  ,  New  York 
(1948),  pp.  239-244. 
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from  and  Pq  by  equating  6  in  (3.  11)  to  -  exp(2f),  (3.  16): 

icp  2ira0  i2irPQ  2ira0  i2Tr(P0  -  j) 

6  =  |fi|e  =  -  e  e  =  e  e 

Hence, 

l  I  a  I2  J/2  2iran  ] 

I6!  =  [ -  - 7TT2]  =  e  and  cp  =  2u(p  -  f)  (3.22) 

(l+A(rT  +  (a(1T 

where  1 6 1  and  cp  are  given  by  eq.  (3.  12).  The  value  of  p^  is  not 
particularly  sensitive  to  the  admittance  function,  but  a.Q  is,  and  in 
particular  is  small  when  the  admittance  is  small.  Note  that  if  A  =  0, 
then  Qq  =  a  as  well.  More  conveniently,  one  can  show  that  the 
admittance  is  related  to  and  Pq  by 

A  =  coth  (3.23) 

In  principle,  four  values  of  |p|  determined  at  four  different 
locations  can  be  used  with  (3.20)  to  determine  the  four  parameters  listed 
above.  If  many  measurements  are  taken,  the  data  can  be  used  in  sets 
of  four,  providing  corresponding  sets  of  the  four  parameters.  The 
results  could  then  be  averaged  over  a  number  of  sets  of  four. 

A  better  way  of  taking  advantage  of  many  data  points  is  to  use  a 
nonlinear  regression  technique  to  fit  the  data  with  the  four-parameter 
representation  (3.20).  A  measure  of  the  difference  between  the  measured 
and  predicted  distributions  is 

E  =  E  [|p(z.)|exp  -  |p(z.;  aQ,  PQ,  K,  P_)|]2  (3.20) 

where  the  sum  is  over  all  data  points.  The  nonlinear  regression 
technique  may  be  used  to  determine  numerically  those  values  of  a^, 

Pq,  K ,  P  which  minimize  E.  This  seems  to  have  been  the  technique 


used  in  Refs.  9-12. 
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Both  of  those  methods  have  been  tried  with  test  problems  defined 
by  specifying  the  four  parameters  first  and  then  using  (3,20)  to  calculate 
"data  points".  Not  surprisingly,  numerical  errors  are  small,  and  the 
parameters  specified  initially  can  be  recovered  with  great  precision. 

The  problem  is  revealed  when  one  follows  the  same  tact,  but 
adds  a  bit  "error"  to  the  specified  data  points.  For  the  sorts  of  values 
one  would  expect  in  the  experiments,  a  1$  error  applied  to  the  data  can 
produce  as  large  as  100$  errors  in  the  values  of  (or  the  real  part 

of  the  admittance  function)  and  *.  Moreover,  it  is  possible  even  to 
extract  values  of  having  the  wrong  sign. 

Several  alternative  procedures  have  been  tried,  including  one  based 
on  the  integral  of  the  amplitude  distribution  and  one  based  on  a  correla¬ 
tion  coefficient.  In  the  course  of  those  calculations,  an  analysis  of  the 
sensitivities  of  the  results  to  errors  in  the  input  data  was  also  carried 
out.  That  served  mainly  to  verify  that  indeed  we  should  have  the  troubles 
we  had  encountered,  although  one  useful  conclusion  was  that  the  inferred 
values  of  the  parameters  should  be  less  uncertain  if  good  data  is  taken 
near  the  minima  of  the  pressure  field.  That  too  is  not  a  startling 
result;  it  follows  simply  from  examining  the  influence  of  the  losses  on 
the  pressure  distribution:  the  influence  is  greater  in  the  vicinity  of  the 
minima  as  shown  by  Figure  20. 

Two  examples  serve  to  show  the  difficulties;  for  simplicity, 

M  =  0.  Figures  20  and  21  show  the  influence  of  changing  by  a 

factor  of  three,  from  -.01  to  -.03  representing  losses  at  the  nozzle. 

The  value  of  K  is  given  by  the  formula  for  viscous  losses  associated 
with  the  classical  acousti  boundary  layer  (Ref.  13,  p.  219) 


} 
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k  =  2/wv  [1  + 


5. 88  x  10 


(3.24) 


where  D  is  the  diameter  of  the  tube  and  the  constant  is  the  value  for 
air  at  one  atmosphere  and  room  temperature.  Note  that  with  all  other 
parameters  fixed,  changing  the  value  of  by  a  substantial  amount 

affects  the  amplitude  and  phase  distributions  only  near  pressure  minima. 

Figures  22  and  23  constitute  a  striking  demonstration  of  the 
crucial  importance  of  measuring  the  phase.  In  this  case,  for  fixed 
values  of  X,  (3q  and  a  small  change  of  amplitude,  changing  the  sign 
of  Oq  causes  essentially  no  discernable  shift  in  the  pressure  distribution 
but  the  phase  is  qualitatively  as  well  as  quantitatively  very  different. 

Even  a  crude  measurement  of  the  phase  distribution  serves  to  distinguish 
between  a  loss  and  gain  of  energy  at  the  end  of  the  tube. 

Finally,  Figures  24  and  25  show  how  insensitive  both  the  amplitude 
and  phase  distributions  are  to  the  value  of  K  except  in  narrow  regions 
immediately  near  the  minima.  It  follows,  as  indeed  we  shall  see  later, 
that  it  is  relatively  difficult  to  obtain  precise  values  for  the  distributed 
losses  from  data  taken  with  an  impedance  tube. 

The  reasons  for  the  awkward  behavior  described  here  can  be 
traced  to  the  form  of  the  amplitude  distribution  (3.20).  Suppose  <  =  0; 
then  the  formula  is 

i~i  *°0r  .2  2  _  ,1/2 

I p  1  =  2P  e  [cosh  ttUq  -  cos  tt(3qJ 

When  data  is  fitted  numerically,  the  amplitude  is  2P  exp(irag),  so 

2 

has  no  influence.  It  then  appears  only  in  cos  ti-Oq  which  is  insensitive 
to  the  sign  of  Qq,  explaining  the  behavior  in  Figure  22. 
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The  sensitivity  may  be  estimated  by  choosing  X  =  0,  {3^  =  0.  5, 
so  |p |  ~  cosh-rra.  Then  the  small  change  of  a  due  to  a  small  change 
of  |p|  is  roughly 

dan  , 

6an  =  — —  6 1  p  | - r~r - 

0  „  )  a  i  1  c  1  it  sinh-rra 

9(Pl 

so  6a  ~  |  p |  /tt  aQ  for  small.  If  a  ~  .01,  Sa^  ~  1 0 1  p  |  and  a 

small  error  in  the  pressure  amplitude  is  amplified  ten  fold  in  its  in¬ 
fluence  on  the  value  of  Oq.  It  is  this  behavior  which  blocks  effective 
use  of  data  for  the  amplitude  distribution  only. 


3.4.2  Numerical  Reduction  of  Data  for  the  Amplitude  and  Phase 
Distributions 

The  technique  we  have  used  to  reduce  the  data  discussed  later  in 
Section  3.  7  is  exactly  that  reported  in  Ref.  18.  In  these  sections  we 
summarize  the  procedure  and  give  the  results  of  several  test  cases 
which  serve  to  verify  the  precision  and  efficiency  of  the  method. 

We  now  take  the  average  flow  to  be  from  left  to  right  towards 
the  nozzle.  Although  we  can  use  the  exact  representation  in  terms  of 
plane  waves,  we  follow  Ref.  18  and  base  the  procedure  on  numerical 
integrations  of  the  linearized  conservation  equations,  written  for  harmonic 


motions , 


— •  du  ,  —  dp  n 
o  t-  +  u  =  0 


i("e  +  P  dl  +  U  d 


(3.25) 


.a  ,  —  du  ,  1  dp  „ a 

iu)u  +  u  -7—  +  —  =  -  Gu 

dz  —  dz 


(3.26) 


-  d& - 2  du  „ 

u  -7s-  +  p  a  — s—  =  0 


(3.27) 


The  parameter  G  is  used  to  introduce  viscous  losses.  It  is  not  dif¬ 
ficult  to  show,  by  substituting  the  earlier  formulas  (3.16)a,b  in  these 
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equations  that  for  G/w  «  1,  G  and  K  are  simply  related: 

G  *  2a  K 

1  +  M 


(3.  28) 


Because  the  density  fluctuation  appears  only  in  (3.25),  the  other 
two  equations  can  be  integrated  as  a  pair.  These  may  be  rearranged 
and  written  in  the  matrix  form, 

^  [Z]  =  [A]  { Z  }  (3.29) 


where 


(3.30) 


[A] 


(iuo  +  G) 
M4 


-  ~  (iit)  +  G) 
u 


(3.31) 


— 2 

Terms  of  order  M  have  been  dropped. 

Now  suppose  that  the  boundary  conditions  are  specified  at 

(zj0  =  (zfoij  = 


z  =  0: 

(3.  32) 


Integration  of  (3.29)  will  lead  to  the  form 

{Z(z)}  =  [T(z)]{Z(0)}  (3.33) 

To  find  the  transfer  matrix  [T],  differentiate  (3.33)  and  compare  with 
(3.29)  to  give 

3“  [T]  =  [A][T]  (3.34) 

with  the  initial  condition  on  [T]: 
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t 


[T(0)]  =  [*  {] 


(3.35) 


The  transfer  matrix  is  then  found  by  integrating  the  pair  of  equations 
represented  by  (3.34).  This  will  then  produce  the  formulas  for  the 
pressure  and  velocity  fields  from  (3.33): 


Pth  =  T,,(z)  p(0)  +  T„(z)  &(0) 


11 


12' 


ath  =  Tzl{z}  P(0)  +  T22(z)  &(0) 


(3.36) 


(3.37) 


The  solution  p^  depends  on  the  parameters  G,  p(0),  u(0).  First 
G  is  fixed  and  p(0)  and  u(0)  are  determined  so  as  to  minimize  the 
complex  quantity  A,  analogous  to  the  error  E  defined  in  connection 
with  the  nonlinear  regression  method: 

4  ■ s  <pth  -  <3-  38> 

The  sum  is  again  over  all  locations  where  data  are  taken.  For  A  to 
be  minimum  with  respect  to  p(0)  and  u(0). 


9A_  -  o 


OA  -  n 


(3.  39) 


dp(0)  "  v  9u(0) 

These  are  two  linear  equations  which  can  be  solved  to  give  p(0),  u(0) 
explicitly  in  terms  of  the  elements  of  the  transfer  matrix  and  the  ex¬ 
perimental  values  for  the  pressure: 

p(0)  =1  [(E  peTn)(E  (T12)2)  -  (E  peT12)(E  t12th)]  (3.40) 


u(0)  =  ^  [EpeT11)(ETuT12)  -  (E  peT12)(E(T11)2)] 


(3.41) 


where 


(3.42) 


-93- 


Note  that  p  ,  standing  for  p  ,  the  measured  pressure,  is  complex. 

©  6Xp 

With  (3.40)  and  (3.41),  the  theoretical  pressure  field  may  now  be  com¬ 
puted  by  using  (3.36). 

The  error  given  by  (3.38)  is  still  a  function  of  the  parameter  G 
which  is  determined  by  iteration,  using  Muller's  method  (Ref.  20)  to 
minimize  |a). 

3.4.3  Some  Test  Cases 

(a)  Sensitivity  of  the  Numerical  Technique 

The  numerical  scheme  described  above  reproduces  specified  exact 

5 

input  to  within  less  than  one  part  in  10  .  That  is,  the  values  of  the 

pressure  field  calculated  for  a  given  admittance  function  at  the  end,  are 

-5 

used  as  "data",  then  the  calculated  admittance  is  within  10  of  the  input 
values  both  real  and  imaginary  parts.  This  conclusion  is  true  over  the 
entire  useful  range  of  values. 

More  importantly,  the  technique  is  relatively  insensitive  to  errors 
or  uncertainties  in  the  data.  This  was  checked  by  specifying  the  viscous 
parameter  G  and  the  admittance  function  and  then  using  as  data  the 
calculated  values  of  the  pressure  field.  Nine  cases  have  been  examined; 
the  results  are  shown  in  Table  10.  The  first  column  lists  the  pairs  of 
Qp,  |3q  used  and  the  second  column  displays  the  corresponding  real 
and  imaginary  parts  of  the  admittance  function.  The  remaining  four 
columns  show  the  computed  values  of  the  admittance  when  the  theoretical 
"data"  are  composed  of  the  exact  values  plus  varying  amounts  of  random 

20.  Muller,  D.A.  "A  Method  for  Solving  Algebraic  Equations  Using 

as  Automatic  Computer",  Mathematical  Tables  and  Other  Aids  to 
Computation  (1956)  pp.  208-215. 
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Table  10 

Results  of  Calculations  for  the  Sensitivity  to  Random  Errors 


Maximum  Value  of  Random  Error 


Boundary 

Conditions 

Theoretical 

Admittance 

±0.  si  P 

max 

±1.  si  p 

max 

±2.  Si  P 

max 

±3.  Si  P 

max 

-0.02004 

-0.02095 

-0. 01933 

-0.  01843 

-0. 01603 

I3I9SBS 

-0.99980 

-0. 10005 

-0. 10054 

-0.99720 

-0. 99483 

-0. 01570 

-0. 01583 

-0.  01537 

-0.  01555 

-0. 01431 

1 

-0. 00018 

-0. 00053 

-0.  00310 

-0.  00172 

-0. 00017 

| 

-0. 04276 

-0.04268 

-0.04236 

-0.04245 

-0. 04175 

1 

+0.99895 

+0.  99817 

+0. 99680 

+0. 99680 

0.  99434 

oo=  0.01 

-0.  05139 

-0. 05240 

-0. 05080 

-0.  04951 

-0. 04701 

Po=  0-25 

-0. 99855 

-0. 99940 

-1.  0040 

-0.99610 

-0. 99326 

a0  =  0.01 

-0.  03140 

-0.03153 

-0.03120 

-0.  03126 

-0.02954 

Po=  0-5 

-0.  00036 

-0. 00072 

-0. 00313 

-0.00193 

-0. 00038 

Qq  =  0.  01 

-0.  07408 

-0.07398 

-0. 07369 

-0. 07355 

-0. 07265 

Po=  0.75 

0. 99711 

0.  99637 

0.  99499 

0.  98920 

0.  99237 

a0=  0.015 

-0. 08233 

-0.08374 

-0. 08223 

-0.  08053 

-0.  07819 

0n  =  0.25 

-0.  99645 

-0.99774 

-1.  0017 

-0.99459 

-0. 99102 

a0=  0.015 

-0.  04708 

-0.04722 

-0. 04697 

-0.04697 

-0. 04488 

0n  =  0.5 

-0.  00054 

-0. 00090 

-0. 00382 

-0. 00213 

-0. 00598 

Qo=  0.015 

-0. 10526 

-0. 10513 

-0.  10489 

-0. 10449 

-0.  10348 

Po?. P^5. 

0.  99430 

0. 99360 

0. 99220 

0. 98662 

0.  98940 

T  T _ _ 1 _ _ _ 1 _ A-  -  —  . _ 1  1 _  _ ^ _  _  1 _  ■ _  a  P  at 


yT 
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errors  distributed  over  the  pressure  distribution.  Even  for  uncertainties 
as  large  as  3.5^  of  the  maximum  pressure,  the  errors  in  the  calculated 
admittance  function  are  less  than  2$.  Thus,  the  amplification  of  un¬ 
certainties  which  accompanies  the  technique  based  on  the  pressure 
amplitude  only,  is  entirely  absent. 

(b)  Results  for  a  Rigid  Flat  Plate 

A  flat  plate  was  installed  as  the  end  condition  and  the  impedance 
tube  operated  without  flow.  In  this  case,  the  admittance  function  is 
ideally  zero  because  the  velocity  normal  to  the  surface  must  vanish. 

Table  11  shows  the  results.  Although  small,  the  admittance  doesn't 
vanish,  a  consequence  of  the  impossibility  of  providing  an  absolutely 
rigid  mounting. 

Also  included  in  Table  1 1  are  the  inferred  values  of  the  attenuation 
constant  <  and  those  calculated  using  the  formula  (3.24).  Note  that 
while  the  admittance  function  varies  only  weakly  with  frequency,  as  it 
should,  the  inferred  values  of  K  do  not  agree  particularly  well  with  the 
predicted  values,  <  ~  J~T.  The  main  reason  for  this  has  been  re¬ 
marked  upon  earlier:  the  distributed  losses  have  mainly  a  relatively 
small  influence  near  the  minima  in  the  pressure  field,  where  the 
experimental  errors  are  greatest.  Hence,  one  should  not  expect  to 
be  able  to  infer  good  values  for  K.  Fortunately,  as  discussion  in 
Section  3.  7  will  clarify  further,  the  values  inferred  for  the  admittance 
function  are  very  insensitive  to  the  value  for  K. 


Admittance 


Frequency  A^  ^exp  (cm  ^  ^th  ^cm  ^  $  Error,  < 


500 

-0. 0056 

0.05902 

0. 0245 

0. 0256 

-  4.  4 1 

700 

-0.0059 

0.  068 

0.  015 

0.0301 

- 100 4 

1050 

-0.  0048 

0.  081 

0.  065 

0. 0381 

+ 100 $ 

(c)  Examples  Showing  Importance  of  Obtaining  Data  Near  Pressure 
Minima 

The  numerical  calculations  discussed  in  Section  3.  4.  1  show  very 
clearly  the  importance  of  obtaining  good  data  near  pressure  minima. 

This  implies  that  best  use  is  made  of  a  given  number  of  data  points  if 
they  are  in  some  sense  concentrated  near  the  minima.  Several  examples 
have  been  calculated  to  show  the  truth  of  this  conclusion,  and  to  demon¬ 
strate  that  one  must  have  data  near  two  minima;  good  definition  of  only 
one  minimum  is  insufficient.  Both  the  end  losses  and  the  distributed 
losses  affect  variations  in  the  stationary  acoustic  field  along  the  axis 
and  good  results  are  obtained  only  if  one  has  some  measure  of  the 
gradients. 

Table  12  shows  the  results  of  one  such  exercise.  A  large  number 
of  data  -  measurements  of  pressure  at  many  axial  locations  -  were  taken. 
The  first  two  rows  show  the  results  for  the  admittance  function  and  the 
attenuation  constant  <  obtained  with  two  different  sets  of  32  points. 
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The  third  row  shows  the  resuits  of  a  calculation  using  only  half  as 
many  data,  but  chosen  from  the  regions  near  two  minima.  The  point 
is  clear,  and  generally  true  for  use  of  the  classical  impedance  tube: 
it  is  best  to  take  data  near  minima,  and  at  least  two  minima  must  be 
well-defined. 

Table  12  An  Example  Showing  the  Importance  of  Data  Near  Pressure 
Minima 


Data  Set 

A<r> 

a<a> 

K 

exp 

*th 

(a)  32  points 

-.0497 

-. 1786 

.  043 

.  038 

(b)  32  points 

-. 0475 

-.  1824 

.  053 

.  038 

(c)  16  points 
near 

minima 

-. 0510 

1833 

.049 

.038 
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3.  5  Summary  of  Tests  Conducted 

Table  13  shows  the  tests  conducted  with  the  two  types  of  nozzles 
illustrated  in  Figure  18.  Because  the  nozzle  entrances  do  not  fair 
smoothly  with  the  lateral  boundary  of  the  impedance  tube,  we  have  no 
possibility  of  comparing  the  results  unambiguously  with  the  well-known 
one-dimensional  theory  given  in  Refs.  6  and  7. 

Each  test  entered  in  Table  12  means  the  acquisition  of  data  for 
both  the  amplitude  and  phase  of  the  pressure  at  various  locations  along 
the  tube,  sometimes  as  many  as  32. 

We  encountered  considerable  problems  with  the  flow  noise  for 
higher  Mach  numbers,  so  most  tests  are  for  M  <  0.  04.  A  total  of 
100  tests  are  entered  in  Table  13.  There  are  of  course  many  not 
included  here  because  they  were  carried  out  to  check  the  apparatus 
and  procedures  for  data  reduction  or  otherwise  were  of  an  exploratory 
nature. 


3 .  6  Calibration  Procedures  and  Estimates  of  Errors 

The  quantities  to  be  determined,  the  real  part  of  admittance  function 
of  the  nozzles,  are  relatively  small  and  their  variation  with  respect  to 
flow  Mach  numbers  and  nozzle  configurations  are  also  relatively  small. 

Thus,  the  measurements  require  considerable  care  and  good  instrumentation 
to  obtain  the  highest  possible  accuracy.  Considerable  effort  was 
devoted  to  achieve  the  necessary  level  of  accuracy. 

3.6.1  Calibration  of  the  Secondary  Air  Supply  System 

A  Brooks  full  view  rotometer  was  used  for  the  calibration 


of  secondary  flow.  As  part  of  this  calibration,  the  rotometer 


was 


Maximum  Minimum 

Submergence  Submergence 


[STS] 


Table  13  Summary  of  Tests  Conducted  for  Measuring  the  Nozzle 
Admittance 


Tests  with  Conical  Nozzles 


f 

.3  kHz 

.  4  kHz 

.5  kHz 

.  6  kHz 

.  7  kHz 

.04 

1 

1 

1 

1 

6 

.  10 

1 

1 

Tests  with  Secondary  Flow  (M 


Submergenc 


2nd  Flow  .  3  kHz 


1 

1 

1 

.  5  kHz 

. 

.  6  kHz 

3 

1 

1 

3 

1 

3 

1 

1 

3 

1 
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separately  calibrated  using  a  Precision  wet  test  meter.  The  flow 
circuit  and  further  details  for  the  calibration  of  the  secondary  flow 
are  given  in  Refs.  3  and  4.  The  supply  pressure  and  flow  rate  are 
controlled  and  maintained  constant  by  a  series  of  coarse  and  fine 
regulators  and  a  fine  metering  needle  valve.  The  pressure  drop  across 
the  porous  plate  in  the  secondary  flow  system  was  measured  with  a 
Statham  differential  pressure  gauge.  The  calibration  of  the  wet  test 
meter  was  provided  by  Precision  Scientific  Company,  while  the  pressure 
gauge  was  calibrated  at  the  Standards  and  Calibration  Laboratory 
of  the  Jet  Propulsion  Laboratory. 

The  results  for  flow  rate  vs.  pressure  drop  across  the  porous 
plate  obtained  in  Refs.  3  and  4  were  adopted  here.  These  porous  plates 
were  chosen  to  provide  very  high  acoustic  impedance,  which  implies 
a  relatively  large  pressure  drop  for  a  given  flow  rate.  The  data 
obtained  are  within  the  tolerances  (±  30$)  quoted  by  the  manufacturer. 

3.  6.  2  Calibration  of  Microphones  and  Other  Instruments 

The  values  of  the  sensitivities  and  the  frequency  response  curves 
for  the  microphones  were  provided  by  the  manufacturer.  Periodic  veri¬ 
fications  were  made  on  the  sensitivity  by  using  a  B&K  piston  phone 
microphone  calibrator.  Model  4220,  used  together  with  a  B&K  measuring 
amplifier.  Model  2607. 

The  Princeton  Applied  Research  lock-in  amplifier  Model  124  A, 
Wavetek  synthesizer  Model  171,  Wavetek  phase  meter  Model  740  and 
the  Data  Precision  digital  voltmeter  Model  3  500  were  sent  periodically 
to  the  calibration  laboratories  of  the  respective  manufacturers  for 
calibration.  The  spectral  Dynamics  tracking  filter  Model  SD-122  and 
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carrier  generator  Model  SD-120  were  calibrated  only  once  at  the  cali¬ 
bration  laboratory  of  the  manufacturer. 

3.  6.  3  Sources  and  Estimates  of  Errors 

Estimates  of  errors  were  carried  out  with  the  best  available 
information  and  some  laboratory  experiments. 

(a)  The  driver  speaker  was  mounted  on  one  side  of  the  tube 
entrance.  Best  results  were  obtained  by  lengthening 

the  impedance  tube  so  that  plane  waves  would  be  achieved 
in  the  test  section.  As  mentioned  earlier,  the  serious 
asymmetries  were  encountered  for  frequencies  above 
750  Hz.,  but  they  have  not  been  accounted  for  in  the 
data  reduction. 

(b)  Determination  of  the  acoustic  pressure  is  affected  by 
the  flexibility  and  mounting  of  the  long  probe.  Spurious 
room  noise  also  has  some  influence.  Neither  of  these 
errors  is  known. 

(c)  The  spacial  resolution  of  the  pressure  measurements  is 
±  0.  1  mm. 

(d)  The  amplitude  linearity  of  the  lock-in  amplifier  enters  in 
determining  the  absolute  values  of  the  signals.  Here,  the 
internal  calibrator  of  the  instrument  was  used  to  estimate 
the  error.  A  practical  difficulty  arises  because  we  need  to 
measure  signals  ranging  from  1#  to  80#  of  the  full  scale 
sensitivity  on  a  single  sensitivity  of  the  lock-in  amplifier. 
When  the  signal  was  only  1#  of  the  full  scale  sensitivity 
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setting,  the  error  in  amplitude  was  determined  to  be 
±  0.  1$,  and  for  signals  greater  than  30$  of  the  full 
scale  it  was  almost  impossible  to  detect  the  very  small 
error. 

(e)  Day  to  day  changes  in  temperature  were  ignored.  They  are 
small  in  the  basement  laboratory  used  for  this  work. 

(f)  The  discharge  coefficient  of  the  nozzles  was  assumed  to 
be  unity.  The  error  in  flow  rate  should  be  less  than  1$. 

3.  7  Results  for  Measurements  of  the  Admittance  Function 

Tests  were  conducted  with  the  two  types  of  nozzles  shown  in 
Figure  18;  the  test  program  covering  100  tests  has  been  summarized  in 
Table  13.  We  shall  discuss  the  results  in  three  groups  for  the  following 
configurations: 

(i)  conical  nozzles  with  no  secondary  flow  (M  =  0.04  and  0.1) 

(ii)  a  nozzle  with  secondary  flow,  minimum  submergence 
(M  =  0.04) 

(iii)  a  nozzle  with  secondary  flow,  maximum  submergence 
(M  =  0.  04) 

3.  7.  1  Conical  Nozzles,  No  Secondary  Flow 

The  inferred  values  for  the  real  and  imaginary  parts  of  the 
admittance  function  are  given  in  Table  14  for  M  =  0.04  and  in 
Table  15  for  M  =  0.1.  Figures  26  and  27  are  plots  of  the  real 
parts.  Repeated  tests  for  M  =  0.04  gave  the  standard  deviations 
listed  and  indicated  by  the  error  bars  drawn  in  Figure  26.  The  ex¬ 
perimental  uncertainty  is  of  the  order  of  10$  except  for  the  series 
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Table  15  Results  of  Tests  of  a  Conical  Nozzle  (M  =  0.  1) 


Frequency 

(Hz) 

A<r> 

A<*> 

2100 

-0.  086 

-0.  96 

1050 

-0. 089 

-0.  3? 

600 

-0.  095 

-0.  187 

300 

-0.084 

-0.072 

Admittance  Function  Measured  for  a  Conical  Nozzle,  M  =  0.  04 


010 


Figure  27.  Admittance  Function  Measured  for  a  Conical  Nozzle 
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at  800  Hz.  ,  for  which  one  of  the  tests  appears  to  be  seriously  in  error, 
for  an  unknown  reason;  given  the  other  results  one  might  guess  that 
that  test  is  likely  invalid. 

If  a  choked  nozzle  behaves  quasi- statically,  the  admittance  function 
is  real  and  given  by  the  formula 

A  =  M  (3.  43) 

n  L 

where  Mn  is  the  Mach  number  of  the  average  flow  entering  the  nozzle. 
This  is  the  ideal  result,  valid  in  the  limit  of  very  low  frequencies. 

For  the  nozzles  used  in  this  program,  the  diameter  of  the  nozzle  entrance 
(see  Figure  18)  is  2.  58  cm,  while  the  diameter  of  the  test  section  is 
5  cm.  Hence,  for  incompressible  flow,  with  Mach  number  M  in  the 
test  section,  M^  =  M  (5.0/2.58)^  =  3.76  M.  The  formula  (3.43)  then 
gives  for  this  short  or  compact  nozzle  approximation, 

0.  150  (M  =  0.  04) 

A  =  _  (3.44) 

0.376  (M  =  0.  1) 

The  values  measured  at  300  Hz.  are  roughly  0.05  for  M  =  0.04  and 
0.085  at  M  a  0.10.  While  the  measurements  show  the  expected  increase 
with  Mach  number,  the  variation  is  less  than  linear  and  the  values  are 
smaller  than  predicted  by  the  simple  theoretical  result. 

The  most  obvious  reason  for  the  discrepancy  is  the  influence  of 
separation  of  the  flow  upstream  of  the  entrance.  This  must  obviously 
happen,  producing  an  effective  nozzle  having  entrance  larger  than  the 
physical  orifice  actually  used. 

We  have  noted  earlier,  we  chose  these  nozzles  because  of  the 
closer  resemblance  to  those  used  for  the  measurements  with  secondary 


flow.  If  a  short  nozzle  is  used  with  entrance  fairing  smoothly  to  the 
test  section,  the  theoretical  values  for  the  admittance  would  be  given 
directly  by  (3.43),  producing  the  values  0.008  for  M  =  0.04  and 
0.02  for  M  =  0.01.  The  measured  results  lie  between  these  and  the 
values  (3.44),  a  reasonable  conclusion. 

That  there  is  considerable  interaction  between  the  acoustic  field 
and  the  regions  of  separation  is  shown  also  by  the  large  measured 
values  for  the  imaginary  parts.  Tables  14  and  15.  There  is  presently 
no  way  to  predict  the  imaginary  parts. 

The  results  therefore  seem  reasonable  so  far  as  we  can  determine, 
and  serve  to  provide  an  estimate  of  the  experimental  errors  one  can 
achieve  with  this  technique.  No  discussion  of  experimental  errors  was 
given  in  Refs.  9-12;  we  therefore  have  no  basis  for  comparison  of  the 
accuracy  of  our  results  with  those  obtained  using  a  blow-down  facility. 

In  general  it  is  better  to  report  the  admittance  as  a  function  of 
the  dimensionless  frequency  s  =  uu  £/a  where  l  is  a  characteristic 
length.  The  usual  choices  for  l  in  the  absence  of  separation  are  the 
length  of  the  convergent  section  or  the  radius  of  the  inlet  section. 

Using  the  latter,  we  have  (r^  =  1.29  cm)  for  the  range  of  dimensionless 
frequency  corresponding  to  the  range  of  real  frequency  300  <  uu  <  2000  Hz. 

0.013  <  s  <  0.  85  (3.  45) 

Because  of  the  smaller  scale  of  our  apparatus,  we  have  been  able  to 
reach  values  of  s  somewhat  less  than  half  those  covered  by  the  work 
in  Refs.  9-12. 

Measurements  of  the  sort  taken  in  this  program  should  ultimately 
provide  a  basis  for  assessing  computations  reported  in  Ref.  21.  Detailed 


calculations  of  the  unsteady  velocity  and  pressure  fields  have  been  performed 
but  the  results  are  presently  unavailable.  In  principle,  it  should  be  pos¬ 
sible  to  compute  the  admittance  function  averaged  over  any  chosen  plane 
upstream  of  the  nozzle,  for  comparison  with  experimental  measurements. 

(See  Section  3.  8  below  for  further  remarks.  )  This  will  be  a  particularly 
important  exercise  for  cases  in  which  the  flow  is  smooth  everywhere.  If 
separat  on  occurs,  as  may  happen  with  secondary  flow,  the  analysis  is 
less  reliable. 

3.7.2  Submerged  Nozzles  with  Secondary  Flow 

Tables  16  and  17  contain  the  values  for  the  admittance  functions 
obtained  with  secondary  flow  and  two  values  of  submergence,  very 
small  but  not  zero,  and  the  maximum  value  possible  with  our  apparatus. 

The  real  parts  of  the  admittance  functions  are  plotted  in  Figures  28 
and  29. 

The  experimental  errors  for  these  tests  increase  greatly  with 
frequency;  the  results  at  2000  Hz.  are  very  doubtful  indeed.  With  the 
caveat  that  there  are  large  uncertainties  attached  to  the  results,  we  may 
offer  the  following  observations: 

(i)  For  frequencies  less  than  1500  Hz.  there  seems  to  be  no 
enormous  influences  of  secondary  flow,  up  to  5^  of  the 
main  flow. 

(ii)  There  also  seems  to  be  no  large  effect  of  submergence,  below 
1500  Hz.  However,  the  maximum  submergence  is  4.3  cm 
and  at  1500  Hz.,  the  wavelength  is  roughly  20  cm,  so  the 
nozzle  spans  less  than  one-quarter  of  a  wavelength. 
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ure  29.  Admittance  Function  Measured  for  a  Nozzle  with  Secondary  Flow, 
Maximum  Submergence,  M  =  0,04 


(iii)  There  is  only  a  vague  suggestion  in  the  data  that  the  real  part 
of  the  admittance  is  increased  by  submergence  when  the  nozzle 
is  submerged  by  more  than  one-quarter  of  a  wavelength.  The 
values  of  the  real  parts  at  2000  Hz.  are,  on  the  average,  less 
for  maximum  submergence. 

The  large  experimental  uncertainties  prevent  us  from  claiming  these 
to  be  credible  conclusions.  they  are  at  best  very  qualitative  and  tentative. 

In  all  cases  for  these  tests,  the  imaginary  part  of  the  admittance 
has  been  determined  to  be  larger  than  the  real  part.  This  conclusion 
is  true  for  all  frequencies,  independently  of  the  submergence  and  value 

of  secondary  flow.  Unfortunately,  the  most  likely  source  of  this 
result  is  not  the  nozzle  or  flow  field,  but  rather  the  system  used  to 
provide  the  secondary  flow. 

As  for  the  flow-turning  measurements  discussed  in  Section  H,  we 
have  been  unable  to  choke  the  porous  plate  surrounding  the  nozzle. 

The  plenum  chamber  is  therefore  not  perfectly  isolated  from  the  acoustic 
field  in  the  test  section  of  the  tunnel.  A  grade  H  sintered  steel  plate 
was  used  to  pass  the  flow.  Its  measured  admittance  has  been  shown 
earlier  in  Figure  9.  For  the  conditions  used  in  these  tests,  the  real 
part  of  its  admittance  is  of  the  order  of  .013,  roughly  1/3  to  1/5  the 
values  measured  for  the  nozzle.  The  imaginary  parts  of  the  admittance 
for  the  plate  have  not  been  determined,  although  there  is  no  reason 
to  believe  that  they  are  larger  than  the  real  parts.  Thus,  the  large 
imaginary  parts  measured  in  the  nozzle  tests  are  almost  certainly  not 
only  a  consequence  of  the  behavior  of  the  nozzle  itself,  but  also  the 
system  used  for  the  gas  supply.  We  cannot  separate  the  one  from  the 
other  with  the  data  we  have  acquired  in  the  present  work. 
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The  major  source  of  difficulty  in  acquiring  data  has  been  flow  noise. 
This  is  of  course  most  troublesome  near  pressure  minima,  exactly  in  the 
regions  where  it  is  most  important  to  obtain  good  data.  The  problem  is 
more  serious  with  the  measurements  of  phase:  that  measurement  carries 
with  it  larger  experimental  error  even  under  the  best  of  conditions. 

Above  Mach  number  0.  1  it  appears  that  the  tracking  filter,  even  with  a 
bandwidth  of  0.  5  Hz.  ,  is  inadequate.  One  would  be  forced  to  use 
machine  processing  of  the  data,  using,  for  example,  correlation  techniques 
or  signal  averaging.  However,  our  experience  suggests  that  with  secondary 
flow  this  becomes  a  very  difficult  problem.  The  noise  we  observed  is  so 
severe  that,  according  to  brief  conversations  with  personnel  at  JPL,  very 
extensive  signal  averaging  would  be  required  to  produce  the  necessary 
effective  bandwidth. 

An  obvious  help  is  to  use  a  stronger  signal,  obtainable  only  with 
pneumatic  drivers.  Since  we  did  not  have  one  at  our  disposal,  we  have 
been  unable  to  verify  the  extent  to  which  this  helps  solve  the  problem. 
Judging  by  the  reports  in  Refs,  9-12,  this  should  help  considerably. 

Just  as  for  the  measurements  of  the  flow  turning  losses,  the 
relative  experimental  errors  can  be  reduced  by  working  with  larger 
scale  apparatus.  The  disadvantage  is  of  course  the  increased  expense, 
arising  from  many  sources.  The  advantages  of  using  as  small  a  facility 
as  possible  are:  low  cost;  ease  of  operation;  and  increased  flexibility. 
However,  if  good  data  cannot  be  obtained,  these  advantages  are  of  little 
value. 

Another  possible  means  of  solving  the  problem  of  data  acquisition 
near  pressure  minima  is  simply  not  to  do  so.  With  the  impressive 
advances  in  laser-Doppler-anemometry  (or  velocimetry)  during  the  past 
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several  years,  it  may  well  be  best  to  return  to  the  more  direct  method 
first  used  in  Ref.  8.  The  admittance  function  can  in  principle  be  deter¬ 
mined  from  data  for  the  pressure  and  velocity  fluctuations  taken  at  a 
position  just  upstream  of  the  nozzle.  This  technique  offers  a  reasonable 
opportunity  for  concentrating  on  acquiring  data  over  one  plane,  thereby 
allowing  one  to  form  the  averaged  quantity  which  is  required  as  the 
boundary  condition  for  stability  analysis.  The  signal/noise  problem  still 
exists,  but  because  the  data  should  be  taken  near  the  nozzle,  it  may  be 
possible  in  some  cases  to  take  data  at  locations  not  near  velocity  or 
pressure  minima.  The  signal  is  greater  and  there  is  no  need  to  define 
pressure  minima. 

3.  8  Remarks  on  the  Plane  of  Definition  for  an  Admittance  Function 

When  data  taken  in  an  impedance  tube  is  reduced  to  obtain  values 
of  the  admittance  function,  the  plane  to  which  the  function  applies  is 
defined  by  the  origin  of  the  coordinate  system  used.  This  is  an  im¬ 
mediate  consequence  of  the  idealized  analysis  described  in  Section  3.  1 
and  of  the  more  general  method  summarized  in  Section  3.4.2.  The 
locations  of  the  pressure  minima  and  maxima;  for  example,  are 
measured  with  respect  to  some  origin  chosen  in  the  experiments. 

The  origin  can  of  course  be  shifted,  and  the  same  data  used  to 
infer  the  admittance  function  for  any  other  plane.  If  there  are  no  losses 
in  the  volume  between  two  transverse  planes,  then  the  admittance  func¬ 
tions  must  be  identical.  We  have  verified  this  property  in  one  simple 
example.  The  major  practical  consequence  is  that  the  plane  chosen  in 
the  experimental  work  must  be  such  as  to  allow  application  of  the 
results  more  generally.  Care  is  required  that  the  plane  selected  should 
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include  downstream  all  influences  which  one  anticipates  finding  in  other 
|  circumstances  to  which  the  results  will  be  applied.  But  if  there  are 

losses  upstream,  they  must  be  accounted  for  in  the  data  reduction. 

For  example,  consider  the  case  of  a  small  conical  nozzle  in  a 
plate,  such  as  that  shown  in  Figure  18.  If  the  origin  (z=0)  is  chosen 
at  the  face  of  the  plate,  influences  of  flow  separation  upstream  should 
not  affect  the  admittance  function.  But  they  certainly  do  affect  the 
acoustic  upstream  where  measurements  will  be  made.  They  must 
therefore  be  accurately  represented  in  the  data  reduction.  In  this 
case  they  cannot  be  modelled  well.  It  is  therefore  necessary  to  place 
the  origin  far  enough  upstream  of  the  face  of  the  plate  that  the  influences 
of  the  separated  region  is  downstream.  The  major  losses  upstream  are 
then  the  viscous  losses  in  the  acoustic  boundary  layer.  The  admittance 
function  inferred  will  then  reflect  the  behavior  of  the  nozzle  and  inter¬ 
actions  between  incident  waves  and  the  separated  regions. 

It  is  especially  important  to  be  aware  of  these  matters  when  dealing 
with  nozzles  exposed  to  secondary  flow.  Because  it  is  not  possible  to 
analyze  such  situations  with  any  confidence,  it  is  essential  that  the  model 
used  for  testing  be  a  faithful  rendition  of  the  entire  aft  end  of  the  motor, 
and  that  the  plane  for  which  the  admittance  function  is  defined  be  carefully 
specified.  The  location  of  that  plane  must  then  be  correctly  incorporated 
in  the  analysis  of  the  stability  of  oscillations  in  the  motor. 
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IV.  SUMMARY  AND  CONCLUSIONS 
4.  1  Flow  Turning  Losses 

The  results  of  this  experimental  program  have  supported  the  main 
qualitative  characteristics  of  the  flow  turning  loss  as  predicted  by  the 
one-dimensional  analysis.  There  is  indeed  a  loss  of  acoustic  energy 
associated  with  average  inward  flow  at  the  boundary,  and  the  loss 
increases  with  the  Mach  number  of  the  average  flow.  There  are  very 
large  uncertainties  associated  with  the  data,  and  the  results  suggest 
that  the  losses  are  several  times  greater  than  those  predicted. 

Owing  to  those  uncertainties  which  reflect  many  difficulties 
arising  in  the  performance  of  the  experiments,  it  is  not  possible  to 
conclude  unequivocally  that  the  actual  losses  are  in  fact  significantly 
greater  than  those  deduced  from  the  approximate  analysis.  Consequently, 
one  should  really  conclude  that  at  the  present  time,  theory  has  predicted 
and  experiment  has  confirmed  the  existence  of  the  losses,  but  their 
quantitative  values  remain  in  doubt. 

This  situation  can  be  improved  in  two  ways;  by  carrying  out  a 
more  elaborate  two-dimensional  analysis;  and  by  performing  further 
experiments.  Theoretical  issues  are  not  our  concern  here;  it  is  cer¬ 
tainly  not  difficult  to  conceive  of  analyses  are  likely  to  improve  upon 
the  one-dimensional  approximation,  although  there  are  no  obvious  simple 
formulations. 

More  pressing  is  the  need  for  better  experimental  work.  It 
appears  that  the  only  way  to  gain  significant  improvement  is  to  use 
larger  apparatus.  The  ultimate  source  of  almost  all  of  our  difficulties 
has  been  the  need  to  infer  a  relatively  small  quantity  as  the  difference 
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between  relatively  large  quantities.  The  flow  turning  loss  per  unit 
area  is  proportional  to  the  Mach  number  of  the  average  flow  and  to 
the  square  of  the  amplitude  of  the  acoustic  field.  Larger  flow  rates 
and  more  intense  acoustic  fields  would  clearly  be  helpful.  And  if  the 
apparatus  is  larger,  the  total  loss  will  increase,  proportional  to  the 
square  of  the  linear  dimension.  With  judicious  design,  it  may  be 
possible  to  achieve  very  substantial  increases  in  the  flow  turning  losses 
relative  to  other  contributions  (such  as  viscous  losses)  thereby  easing 
the  problems  of  measurement. 

4.  2  Measurements  of  the  Admittance  Function  for  a  Choked  Nozzle 

The  unfortunate  problems  with  data  processing,  discussed  in 
Section  3.4,  eventually  prevented  us  from  gathering  satisfactory  data 
for  the  admittance  function.  The  main  achievements  have  been  to 
verify  that  the  method  developed  at  Georgia  Tech  is  a  thoroughly 
satisfactory  method  for  reducing  data  taken  in  an  impedance  tube  with 
flow;  and  the  technique  based  on  using  a  continuously  operating  super¬ 
sonic  wind  tunnel  will  work,  but  important  improvements  must  be  made. 
There  is  no  need  to  discuss  data  processing  further. 

Owing  to  the  limited  data  eventually  taken  in  this  program,  no 
unusual  or  new  general  conclusions  can  be  drawn.  The  experimental 
errors  are  not  excessively  large  for  the  cases  run  with  no  secondary 
flow:  roughly  20#  or  less  for  the  real  part  of  the  admittance  function 

of  a  small  conical  nozzle  in  a  plate.  With  secondary  flow,  the  errors 
were  unacceptably  large,  in  some  cases  larger  than  the  average  values. 
The  errors  increase  both  with  flow  rate  and  with  frequency,  for  reasons 
which  we  do  not  comprehend  at  this  time.  Because  of  the  large  errors, 
it  is  not  possible  to  offer  any  general  conclusions  concerning  the 
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influences  of  secondary  flow  and  submergence.  The  data  show  only 
that  the  effects  are  not  accountable  as  factors  of  ten. 

By  far  the  major  experimental  problem  is  signal/noise  at  pressure 
minima.  Clearly  the  sensible  solution  is  to  use  a  larger  source  for  the 
acoustic  field.  Electromechanical  speakers  are  probably  not  satisfactory, 
but  there  are  at  least  three  other  possibilities:  pneumatic  drivers  or 
rotary  valves;  crystal  drivers  (which  seem  to  be  difficult  to  find  for  the 
audio  range);  and  oscillating  obstructions  in  the  flow.  The  last  is 
really  a  form  of  rotary  valve  in  the  broad  :st  sense.  Increasingl  the 
signal  is  a  soluble  problem,  although  because  the  noise  increases 
rapidly  with  increasing  Mach  number,  only  exploratory  tests  will 
verify  the  limits  to  which  one  may  go. 

Use  of  a  continuously  movable  microphone  to  take  data  is  essential, 
but  in  future  work,  an  arrangement  should  be  designed  to  make  the 
measurements  in  the  lateral  boundary.  Although  we  minimized  the 
interference  caused  by  a  microphone  concentrically  placed  in  the  tube, 
there  is  clearly  some  influence.  Much  worse,  the  necessity  to  take 
data  through  a  long  tube  greatly  reduces  the  total  signal  and  aggravates 
the  signal/noise  problem. 

With  a  continuously  operating  system,  it  is  relatively  straight¬ 
forward  to  study  the  effects  of  secondary  flow.  It  is  possible  not  only 
to  determine  the  influences  on  the  acoustical  properties,  but  the  possibility 
exists  for  making  detailed  observations  of  the  flow  field.  That  is  an 
important  subject  which  we  did  not  pursue  in  the  present  program  beyond 
a  small  number  of  exploratory  tests.  The  basic  problem  of  the  structure 
of  the  acoustic  field  in  an  average  flow  field  approaching  fully  developed 
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pipe  flow  remains  untouched.  We  have  obtained  no  measurements  of 
the  flow  field  in  the  vicinity  of  a  nozzle  with  secondary  flow,  a  subject 
which  merits  considerable  attention.  It  appears  that  this  problem  would 
benefit  greatly  from  application  of  laser-Doppler  velocimetry.  In 
principle  one  should  be  able  to  obtain  measurements  of  both  the  fluctuat¬ 
ing  pressure  and  velocity  as  functions  of  position  in  regions  upstream 
of  the  nozzle.  The  data  would  provide  a  means  of  determining  the 
admittance  function  directly,  both  as  a  function  of  position,  and  averaged 
over  a  plane.  This  would  be  a  helpful  independent  check  of  the  results 
obtained  using  the  apparatus  as  an  impedance  tube.  The  device  merits 
development  as  a  means  of  making,  eventually,  routine  measurements 
of  the  behavior  of  choked  nozzles. 


V.  REFERENCES 


Culick,  F.  E.  C.  "Remarks  on  Entropy  Production  in  the  One- 
Dimensional  Approximation  to  Unsteady  Flow  in  Combustion 
Chambers",  Comb.  Sci.  and  Tech.  .  V.  15  (197?)  pp.  93-97. 

Culick,  F.  E.  C.  "The  Stability  of  One-Dimensional  Motions  in  a 
Rocket  Motor",  Comb,  and  Sci.  and  Tech.,  V.  7  (1973)  pp. 

165-175.  “  ~  "  . . *  ’ 

Magiawala,  K.  "Measurements  of  Energy  Exchange  Between 
Acoustic  Fields  and  Non-Uniform  Flow  Fields",  Ph.  D.  Thesis, 
California  Institute  of  Technology  (May  1978). 

Culick,  F.  E.  C.  and  Magiawala,  K.  ,  "Measurements  of  Energy 
Losses  Associated  with  Interactions  Between  Acoustic  Waves 
and  a  Steady  Flow  Field",  Air  Force  Rocket  Propulsion  Laboratory, 
Report  AFRPL-TR- 78-6  (March  1978). 

Magiawala,  K.  and  Culick,  F.  E.  C.  "Measurements  of  Energy 
Exchange  Between  Acoustic  Fields  and  Non-Uniform  Steady 
Flow  Fields",  J.  Sound  and  Vib.  (to  be  published). 

Tsien,  H.  S.  "The  Transfer  Function  of  Rocket  Nozzles",  ARSJ, 

V.  22,  No.  3  1  (May- June  1952),  p.  139. 

Crocco,  L.  and  Cheng,  S.  -T.  Theory  of  Combustion  Instability  in 
Liquid  Propellant  Rocket  Motors,  AGARDcgraph  No.  8,  Butterworths 
Publications,  Ltd.,  London  (1958)  Appendix  C. 

Crocco,  L.  ,  Monti,  R.  and  Grey,  J.  "Verification  of  Nozzle  Ad¬ 
mittance  Theory  by  Direct  Measurement  of  the  Admittance 
Parameter",  ARSJ,  V.  31,  No.  6  (June  1961)  pp.  771-775. 

Zinn,  B.  T.  et  al  "Experimental  Determination  of  Three-Dimensional 
Liquid  Rocket  Nozzle  Admittances",  AIAA  J.  ,  V.  11,  No.  31 
(Mar.  1973),  pp.  267-272. 

Bell,  W.  A.,  Daniel,  B.  R.  and  Zinn,  B.  T.  "Experimental  and 
Theoretical  Determination  of  the  Admittances  of  a  Family  of 
Nozzles  Subjected  to  Axial  Instabilities",  J.  Sound  and  Vib.  , 

V.  30,  No.  2  (1973),  pp.  179-190. 

Zinn,  B.  T.  et  al  "Damping  of  Axial  Instabilities  by  the  Minuteman  II, 
Stage  III  and  Minuteman  III,  Stage  III,  Exhaust  Nozzles",  Georgia 
Institute  of  Technology,  Report  AFRPL-TR- 72-71  (August  1972). 

_  "Nozzle  Design  Considerations  for  Attenuation  of 

Axial  Instabilities  in  the  Minuteman  II  and  III,  Stage  III  Rocket 
Motors",  Georgia  Institute  of  Technology,  Report  AFRPL-TR-73-69 
(September  1973). 


j 


-123- 


13.  Culick,  F.  E.  C.  (Ed.),  "T-Burner  Testing  of  Metallized  Propellants", 
AFRPL  Report  TR- 74-28  (Oct.  1974). 

14.  Shapiro,  A.  H.  The  Dynamics  and  Thermodynamics  of  Compressible 
Fluid  Flow.  The  Ronald  Press  Company,  New  York  (1953). 

15.  Rauscher,  M.  Introduction  to  Aeronautical  Dynamics.  John  Wiley 
and  Sons,  Inc.,  New  York  (1953)  p.  50. 

16.  Van  Moorhem,  W.  K.  "An  Investigation  of  the  Origin  of  the  Flow- 
Turning  Effect  in  Combustion  Instability",  17th  JANNAF  Combustion 
Meeting  (September  1980). 

17.  Culick,  F.  E.  C.  "Linear  Analysis  of  One-Dimensional  Oscillations 
in  a  Variable  Area  T-Burner",  9th  JANNAF  Combustion  Meeting 
(1972). 

18.  Baum,  J.  D.  ,  Daniel,  B.  R.  and  Zinn,  B.  T.  "Determination  of  Solid 
Propellant  Admittances  by  the  Impedance  Tube  Method",  AIAA  18th 
Aerospace  Sciences  Meeting  (Jan.  1980)  AIAA  Paper  80-0281. 

19.  Morse,  P.  Vibration  and  Sound,  McGraw-Hill  Book  Co.  ,  New  York 
(1948),  pp.  239-244. 

20.  Muller,  D.  A.  "A  Method  for  Solving  Algebraic  Equations  Using 
as  Automatic  Computer",  Mathematical  Tables  and  Other  Aids  to 
Computation  (1956)  pp.  208-215. 


